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ABSTRACT 
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This  thesis  describes  the  development  of  a  catalogue  of 

waverider  shapes  generated  from  axisymmetric  conical  flow. 

A  simple  but  reliable  method  of  calculating  the  viscous  drag 

is  presented  and  over  one- hundred  configurations  are  shown 

with  various  aerodynamic  and  pertinent  geometric  design 

factors.  All  of  the  configurations  presented  are  based  on 

Mach  number  of  l1*)  and  Reynolds  number  (based  on  the  length  of 

—  jo  L>i‘-  vv<  # 

the  waverider)  of  .10^  However,  t he*~ana  1  ys i s  presented  here 
adapts  itself  to  any  Mach  or  Reynolds  number.  The  lift-to- 
drag  ratios  presented  range  from  4.32  to  5.06  while  the 
average  skin  friction  coefficient  ranged  from  0.00176  to 
0.00205.  The  volume  to  area  ra tid-^2/3/Sp  Ranged  from  0.224  to 
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Nomenclature 


VAn 


(L/D), 


ratio  of  base  areas  of  an  arbitrary  wave- 
rider  with  the  idealized  cone  derived 
waver ider 

aspect  ratio,  base  of  waverider  divided 
by  the  length  of  waverider 
drag  coefficient 

average  skin  friction  coefficient  for  the 
compression  surface 

skin  friction  coefficient  based  on  conditions 
just  outside  the  boundary  layer 
average  skin  friction  coefficient  for  the 
freestream  surface 

average  skin  friction  coefficient  for  the 
compression  surface 

average  total  skin  friction  coefficient 
same  as 

thermal  conductivity  at  the  wall 
length  of  basic  cone 
length  of  waverider 
viscous  lift  to  drag  ratio 

Mach  number  just  outside  the  boundary  layer 
Mach  number  at  infinity  or  freestream 


v  i  i  i 


condi tions 


Pr  Prandtl  number 

q  dynamic  pressure 

qw  heat  flux  at  the  wall 

r  distance  from  the  origin  of  axis  to  a  point 

on  the  waverider  surface 

Rcb  distance  from  origin  to  an  arbitrary  point  on 

the  the  compression  surface  projected  to  the 
base  plane 

R  e  Reynolds  number  based  on  conditions  just  out-side 

the  boundary  layer 

Rew  Reynolds  number  based  on  freestream  conditions 

R^  distance  from  origin  to  an  arbitrary  point  on 

the  freestream  surface  projected  to  the  base 
p  1  ane 

Sp/Spi  ratio  of  planform  areas  of  waverider  to  the 

idealized  waverider 

S/S  ratio  of  wetted  area  to  planform  area 

SSD  shock  standoff  distance  in  the  base  plane 

Sj  Stanton  number 

wetted  area  of  compression  surface 

wetted  area  of  freestream  surface 

adiabatic  wall  temperature 

temperature  at  the  wall 


volume  ratio  of  waverider  to  idealized  wave- 
nond imens iona 1 i z ed  volume  to  planform  area 
ratio 

halt  angle  of  shock  cone 
halt  angle  ot  basic  cone 
boundary  layer  thickness 

boundary  layer  thickness  in  the  transformed 

coord i na  tes 

angle  between  -<£]  and 

dihedral  angle  ot  waverider  in  the  base  plane 
dihedral  angle  ot  waverider  at  an  arbitrary 
c rossec  tion 

ratio  ot  specific  heats 

viscosity  coefficient  just  outside  the 
boundary  layer 

viscosity  coefficient  at  the  wall 
3/6,  radius  of  shock  cone  in  the  base  plane 
of  nondimensiona 1  coordinate  system 
radius  of  shock  cone  at  an  arbitrary  cross 
section  of  nond imens l onal  coordinate  system 
surface  stresses  at  the  wall 


CHAPTER  1 


Introduction 

1.1  Introductory  Remarks 

Lifting-body  configurations  are  currently  of  interest  as 
a  means  for  achieving  long-range  and  high-performance  hyper¬ 
sonic  missile  and  aircraft  characteristics.  Semi-empirical 
methods  used  heretofore  for  calculating  the  flow  fields  and 
aerodynamics  of  blended  wing-body  shapes  are  not  appropriate 
for  modern  hypersonic  flight.  Design  requirements  for 

I  " 

hypersonic  aircraft  are  discussed  by  Townend  and  Stollery*. 
Elaborate  computer  codes  have  become  necessary  for  predicting 
the  charac ter istics  of  arbitrary  configurations.  Since  these 
codes  are  still  in  evolutionary  stages  and  are  quite 
expensive  as  well  as  complex,  a  need  for  a  simplified 
calculation  scheme  exists. 

Waverider  analysis  provides  one  means  for  studying 
blended  body  configurations  on  a  relatively  simplified 
scale.  The  work  presented  here  is  well  adapted  to  the 
personal  computer.  Shapes  derived  by  this  analysis  yield 
flow  fields  and  aerodynamic  properties  that  are  known  for 
their  particular  on-design  conditions,  and  allow  for 
systematic  variation  of  parameters  pertinent  to  the  design 
process . 
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An  historical  background  on  waveriders  is  provided  by 
Kuchemann^  and  Hemsch  and  Nielsen'  offer  additional 
perspectives.  Also,  in  recent  times  there  have  been  several 
papers  dealing  with  geometry  variations,  optimization, 

viscous  effects,  experimental  verification,  and 

5-15 

appl ications . 

1.2  Primary  Goals 

This  study  has  two  primary  goals.  The  first  is 
essentially  an  extension  of  the  work  by  Rasmussen  and 
Stevens".  In  this  sense,  this  work  is  a  catalog  of  many 
shapes  and  variations  of  waveriders  (over  100  are  presented) 
that  can  be  generated  by  axisymmetric  flow  past  a  cone. 
Lift,  drag,  and  friction  coefficients  as  well  as  lift-to- 

drag  ratios  are  C3lculated  together  with  pertinent  design 

factors  such  as  projected  planform  area,  base  area,  wetted 
surface  area,  volume  ratio,  and  aspect  ratio. 

The  other  goal  of  this  study  is  the  calculation  of  the 
viscous  drag  in  a  simple  but  reliable  fashion.  The  viscous 
drag  has  been  calculated  previously  by  an  average  skin 
friciton  coefficient,  such  as  outlined  by  reference  8,  or  by 
means  of  an  elaborate  integral  boundary- 1  ayer  method1^.  Here 
however,  a  new  method  is  presented  for  laminar  boundary 

layers  that  allows  for  variations  in  Mach  number,  Reynolds 


3 

numbers  based  on  body  length,  and  wall  temperature 
conditions.  The  method  involves  simple  quadratures 
associated  with  the  wetted  surface  area  of  the  waveriders. 

It  is  hoped  that  we  can  present  some  configurations  that 
heretofore  were  not  envisaged  using  cone  derived  waveriders. 
None- the- less  these  perhaps  unusual  shapes  still  render  flow 
fields  that  are  known  based  on  supersonic  flow  past  a  cone. 


CHAPTER  2 


Development  of  Waverider  Shape 

2.1  Introductory  Remarks 

We  will  start  by  defining  a  coordinate  system  and  basic 
parameters.  Then  we  will  determine  in  a  systematic  way  the 
waverider  shape  starting  with  the  base  plane  view.  Finally, 
we  will  present  views  at  arbitrary  cross  sections  and  a  view 
from  overhead. 

2.2  Coordinate  System  and  Basic  Definitions 

To  begin  we  will  use  the  spherical  coordinate  system  shown 
below  in  figure  2-1. 


Z 


Figure  2-1  Spherical  Coordinate  System 
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This  figure  lends  itself  to  the  following  set  of  basic 
equations . 

x  =  rsin0cos$,  y  *  rsin@sin$,  z  =  rcos©  (2-1) 
We  now  wish  to  impose  that  9  be  restricted  to  small  angles 
and  we  reduce  this  set  of  equations  to, 

x  =  r@cos#,  y  =  r©sin<|>,  z  =  r  (2-2) 
Next,  we  rotate  the  conventional  spherical  coordinate  system 
and  impose  a  cone  symmetric  about  the  z-axis  with  the  vertex 
at  the  origin,  and  establish  a  free  stream  Mach  number  aligned 
with  and  pointing  in  the  +  z  direction.  This  creates  a  shock 
cone  that  is  also  symmetric  about  the  z-axis.  Figure  2-2 
below  shows  this. 
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We  now  want  to  construct  a  non— conical  body  that  will  have 
a  flow  field  that  remains  conical  between  the  body  compression 
surface  and  the  shock  surface.  This  is  done  using  the 
procedure  in  reference  14. 

To  begin,  we  define  the  trailing  edge  of  the  free  stream 
surface.  In  this  thesis  we  will  use  a  sixth-order  polynomial 
to  describe  this  line.  Now,  we  trace  undisturbed  free 
streamlines  (i.e.,  streamlines  parallel  to  M,)  upstream  until 
we  intersect  the  shock  cone,  thus  describing  the  freestream 
surface.  Next,  we  trace  corresponding  streamlines  through  the 
shock  layer  until  we  intersect  the  base  plane  of  the  cone, 
thus  describing  the  compression  surface  of  the  waverider. 

When  we  determine  the  streamlines  in  the  compression 
surface  that  correspond  to  those  in  the  freestream  surface  we 
require  the  aid  of  figure  2-3.  Imagine  that  we  pass  a  plane 
through  the  axis  of  the  basic  cone.  This  plane  then 
identifies  cor responding  streamlines  for  the  freestream  and 
compression  surfaces.  For  example,  we  may  choose  any  point 
on  the  trailing  edge  of  the  freestream  surface  and  pass  a 
plane  through  the  axis  of  the  basic  cone  which  also  intersects 
this  point.  Now  staying  in  the  plane  we  trace  upstream  until 
we  intersect  the  shock.  Then  we  trace  back  a  streamline  that 
intersects  the  shock  and  stays  in  the  plane  until  we  get  back 
to  the  base  of  the  cone.  Now  if  we  sweep  the  plane  from  point 
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to  point 

along  the 

f reestream 

surface 

we  define 

the 

compression 

surface . 

Theref ore , 

simply  by 

describing 

the 

trailing  edge  o  f  the  freestream  surface  we  have  means  to 
determine  the  entire  shape  of  the  waverider  body. 
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According  to  reference  14,  the  freestream  and  compression 
surfaces  are  described  by 

r$  »  rf($)|3  (2-3) 

r(©c2  -  S2)l/2  »=  rs($)((32  -  S2)l/2  (2-4) 

where  r  is  a  line  from  the  origin  of  our  coordinate  system 
to  any  arbitrary  point  on  the  freestream  or  compression 
surface,  ©  is  the  angle  between  r  and  the  axis  of  the  cone 
(z-axis),  and  r?  is  the  line  from  the  origin  to  the 
intersection  of  the  compression  and  freestream  surfaces. 
Also,  (3  and  6  are  the  half  angles  for  the  shock  cone  and 
the  basic  cone  respectively.  Figure  2-4  below  shows  a  2-D 
side  view  of  these  definitions. 


Figure  2-4  Two  Dimensional  Side  View 
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2.3  The  base  plane  view 

We  would  like  to  represent  the  trailing  edge  of  the  free- 
stream  surface  by  a  4-term,  sixth-order  polynomial  like  the 
foil  owing 

x  =  A  +  By2  +  Cy4  +  Dy6  (2-5) 
A  representation  of  this  form  will  allow  us  to  be  symmetric 
about  a  plane  passing  through  the  z-axis  in  the  x-direction. 
Consider  a  base  plane  representation  as  is  shown  below  in 
figure  2-5  along  with  the  above  equation  (2-5). 


Figure  2-5  Baseplane  Representation  (dimensional) 

We  can  see  from  equation  (2-5)  that, 

x  =  A  when  y  *  0  (2-6) 

Also,  considering  figure  2-4  we  have  the  distance  from  the 
origin  to  the  shock  as  lsinft  .  Note  that  this  is  also  the 


distance  from  the  origin  to  the  intersection  of 


the  two 
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surfaces  since  the  shock  remains  conical.  This  perhaps 
suggest  that  we  normalize  with  lsinfl  .  However,  considering 
the  relationship  between  fl  and  6  given  in  reference  14  as, 

* „  jl  .  [(mi  t  _i_  ,r 

a  L  2  (».«n 

we  choose  to  normalize  with  IsinS  or  16  since  we  have 
restricted  ourselves  to  small  angles.  Dividing  equation  (2- 
5)  by  16  yields 

x /  (  1  6  )  =  A/ (16)  +  1/(16  )CBy2  +  Cy4  +  Dy6]  (2-8) 

Now  as  a  matter  of  convenience  since  8,  C,  and  D  are  as  yet 
arbitrary  we  will  define  them  in  the  following  way 

B  =  b2/(16),  C  »  b4/(16)3,  D  =  b6/(16)5  (2-9) 

where  bj,  b^,  and  bfc  are  yet  to  be  described.  Also,  we  will 
define  A/ (16)  as  a  new  parameter  with  the  symbol  RQ  . 
Doing  this  we  rewrite  equation  (2-8)  as, 

X  =  RQ  +  b2YZ  +  b4Y4  +  b^Y6  (2-10) 

where , 

X  =  x/  ( 1 6  )  ,  Y  »  y/ (  1 8 )  ( 2-11 ) 

This  normalized  system  renders  the  following  new  baseplane 
representation  shown  in  figure  2-6.  It  is  this  new  normalized 
system  that  we  will  use  to  develop  a  method  of  determining  an 
infinite  variety  of  waverider  baseplane  shapes. 

This  system  allows  the  following  conditions: 


i  )  X  =  Rq  when  V  =  0 
ii)  X  =  ccos$)  when  Y  *  asinfj 
Additionally,  we  will  require  that  the  slope  at  the  tip  of  the 
waverider  be  between  zero  and  This  helps  to  prevent 

double  intersections  with  the  compression  surface.  We  will 
call  this  condition  three, 
iii)  0  <  dX/dY  < 

The  four  terms  of  equation  (2-11)  are  determined  using  these 
conditions.  Substitute  condition  ii)  into  (2-11)  to  yield, 
tfcos$|  *  Rp+bjO^sin2# j+b^c^sin^j+b6a6sin6f j  (2-12) 


Figure  2-6  Baseplane  Represen ta tion  (non-dimensional) 
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We  see  that  in  order  for  (2-12)  to  hold  we  must  choose 
R  ,  $j,  and  b^  accordingly.  To  do  this  we  may  choose 
any  one  variable  to  be  dependent  and  let  the  others  remain 
independent . 

For  example,  choose  bj  to  be  the  dependent  variable. 
Solving  equation  (2-12)  for  b2  yields, 


crcoslij-Rg-b^c'sin^j-bjCr^sin4^] 

5  5 

a  sin 


(2-13) 


We  may  therefore  use  (2-10)  to  describe  a  general  waverider 
freestream  edge  providing  we  use  (2-13)  to  determine  bj  once 
we  have  arbitrarily  chosen  Rq,  $j,  bp  and  bft. 

Additionally,  we  will  restrict  attention  to  the  cases  when 
condition  (iii)  holds  as  well.  To  see  if  (iii)  holds  we 
merely  calculate  dX/dY  =  2bjY  +  4b^Y^  +  6b^  and  see  if  the 
criterion  is  met. 

We  now  have  a  means  to  make  various  freestream  trailing 
edge  shapes.  Referring  to  the  above  coordinate  system  we  see 
that  the  negative  most  value  of  Y  is, 

%  =  ~asin*l  (2-14) 
and  clearly  the  positive  most  value  of  Y  is, 

Y^  =  -t-csin^j  (2-15) 
All  other  values  of  Y  must  lie  in  between  (2-14)  and  (2-15). 

We  therefore  start  at  Y  =  -asin$j  and  march  away  in 


13 


small  increments  until  we  reach  Y  =  ♦csin^  .  At  each  step 
along  the  way  we  stop  and  use  (2-10)  to  determine  X.  We  now 
have  many  X,  Y  data  points  available  to  plot  the  freestream 
edge.  If  the  increments  are  significantly  small  we  may  use 
a  straight  line  approximation  in  between  data  points. 

The  above  is  a  procedure  to  determine  the  freestream  edge 
in  the  base  plane.  We  now  develop  the  procedure  to  determine 
the  compression  surface.  Recall  figure  2-6.  At  each  set  of 
points  X,  Y  we  may  calculate  R  .  by  using  the  Pytnagorean 
t  heorem 

R.b  =  (X2  +  V2)1/2  (2-16) 
where  R>()  is  the  distance  from  the  origin  of  the  base  plane 
coordinate  system  too  a  point  on  the  surface  of  the  freestream 
edge  in  the  base  plane. 

Now  wo?  must  determine  the  compression  surface  that  goes 
with  a  chosen  freestream  surface.  This  is  accomplished  in 
reference  14  and  is  given  as 

Rcb  “  C1  *  (  ff2  "  DR2^/^2]172  (2-17) 

We  will  now  define  the  X,  Y  data  points  for  the 
compression  surface  as  X£  and  Y^  to  prevent  confusion.  To 
determine  the  Xc,  Y£  data  points  for  the  compression  surface 
we  use  simple  trigonometric  relations, 

Xf  =  Rc(jcos$,  Yf  =  R^sin#  (2-18) 
But,  cos#  =  X / R^  and  sin#  =  Y/R^  and  we  rewrite  (2-18)  as 
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xc  =  *<Rcb/FV’  yc  =  v(Rcb/R.b>  (2~19) 

The  above  gives  a  systematic  way  to  explore  virtually  an 

infinite  number  of  waverider  base  plane  shapes.  All  we  need 

do  is  choose  different  variables  subject  to  the  given 

conditions  and  plot  our  picture. 

We  will  now  look  at  arbitrary  crossections  of  the 

waverider  viewed  from  a  base  plane  perspective. 

2.4  Arbitrary  Crossectional  Views 

The  freestream  surface  at  any  arbitrary  cross  section  is 
still  given  by  the  general  equation, 

X  =  R0  +  b2Y2  *  b/  +  b/  (2-20) 

however  the  range  of  V  is  clearly  not  the  same.  Recall  figure 
2-3.  From  the  figure  we  can  see  that  <$z  is  dependent  on 
location  along  the  Z  axis.  Note  is  the  angle  between  the 
Z-axis  and  the  outermost  value  of  Y. 

Now  recall  condition  ii)  given  earlier.  We  can  use  the 
same  condition  at  any  arbitrary  cross  section  simply  by 
replacing  with  4>z  and  o  with  a z.  Thi'a  gives, 

X  =  (jjCOS^j  when  Y  =  OjSin^  (2-21) 

substituting  (2-21)  into  (2-20)  gives, 

cr2cos$z  =  R^bjCr^sin^+bjCT^sin^+b^a^sin6^  (2-22) 

There  are  two  unknowns  in  (2-22).  They  are  cz  and  4>z . 
Considering  figure  (2-7)  below  we  can  produce  a  relationship 
that  will  give  us  a  . 


Figure  2-7  Side  View  (non-dimensional) 

Recall  that  <rz  is  the  distance  from  the  Z  axis  to  the  line 
of  intersection  of  the  shock  and  waverider  body  (in  the  non- 
dimensional  system).  But,  this  is  also  the  distance  from  the 
Z  axis  to  the  shock  in  the  X-direction.  Therefore  simple 
triangle  relations  give  us, 

at  =  R0  +  Z/l^o  -  Rq)  (2-23) 
where  Z^/l^  is  the  percentage  back  from  the  leading  edge, 
i.e.,  0.1,  0.35  etc... 

Now  equation  (2-23)  gives  us  az,  and  then  (2-22)  will  give 
us  Due  to  the  difficulty  of  solving  (2-22)  for  we  may 
use  a  numerical,  instead  of  analytical  approach  to  solving 
it.  Once  and  az  are  known  the  free  stream  edge  is  simply 
plotted  as  before  by  replacing  a  and  by  a2  and  $2  in 


\ 
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equations  (2-14)  and  (2-15).  We  may  call  these  equations  (2- 
14)  and  (2-15)  modified. 

We  now  proceed  to  determine  the  compression  surface  at  an 
arbitrary  cross  section.  Recall  equation  (2-4)  describing 
the  compression  surface, 

r2i0c2  -  62)  =  r2s($)((32  -  $2)  (2-24) 

Dividing  through  by  1  we  obtain, 

(r2/l2)<ec2  -  S2)  =  r2s($)/l2(02  -  62)  (2-25) 

Now  we  define  as  the  angle  between  r  and  the  Z-.  xis 

where  r  is  a  line  from  the  origin  to  any  point  on  the  free- 
stream  trailing  edge  in  the  base  plane.  With  this  definition 
then  we  can  say  that, 

10*  =  r5p  (2-26) 

where  r  s  is  a  line  from  the  origin  to  the  point  of 
intersection  of  the  two  surfaces.  Dividing  through  by  6  we 

obtain , 

10*/  6  =  rf3/S  =  rso  (2-27) 

and  solving  for  0*/6  yields, 

©*/€  =  rtc/l  (2-28) 

Now  0*/S  is  also  equal  to  R*  and  we  can  write, 

R*  =  ria/1  (2-29) 

or 

rs2/l2  =  R*2/j2  (2-30) 

Substituting  (2-30)  into  (2-25)  gives, 


(  r2/  1 2)  ( ©c2  -  S2)  =  (R^/o2)  ( |32  -  62)  (  2-31  ) 

2 

and  dividing  through  by  S  yields, 

r2/l2(@c2/62  -  1)  *  R#b2/a2(  (32/  62  -  1)  (2  -  32) 

Solving  (2-32)  -for  ©C2/S2  and  replacing  r  with  z  yields, 

i 

ec2/52  =  1  +  (l2/z2)[(c2  -1  (2-33)  j 

We  note  that  ©C2/S2  =  R£l)2  and  therefore  (2-33)  is  remarkably 

2  2 

close  to  equation  (2—17)  except  for  the  factor  1  /z  » 

i 

Since  we  want  to  view  the  arbitrary  cross  section  from  a 
base  plane  perspective  we  want  to  project  ©c  onto  the  base 
plane.  Figure  2-8  below  represents  this. 
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We  see  that, 

z  9  =  1  (  @  ) 
c  1  c'proj 

(2-34) 

or 

(  9r  )  =  29  /  1 

c  proj  c 

(2-35) 

Substituting  (2-35)  into  (2-33)  yields, 

=  z2/  1 2  +  (cr2  “I  iR^t* ) /a2 

(2-36) 

From  the  above  figure 

we  see  that, 

2.  *  2  "  Ro1/a 

(2-37) 

and 

lS  " 

1  -  R0l/a  -  1(1  -  Rq/<t) 

(2-38) 

Equation  (2-38)  lets 

us  normalize  (2-36)  with  lw. 

Solving 

(2-38)  for  1  and  substituting  into  (2-37)  we  get  after 
rearranging , 

zu/l9  -  lz/1  -  RQ/a)/(l  -  R/ff)  (2-39) 

We  now  solve  (2-39)  for  z/\  to  yield, 

z/1  -  (  z/l„)  ( 1-  R0/ct)  RQ/ff  (2-40) 

Now  substitute  (2-40)  into  (2-36)  to  give  the  desired  result, 
(Rd,2)proj  *  (z/lMU-R0/o)+R0/a)2+(a2-l)R#b2(#)/cr2  (2-41) 

We  now  determine  X£,  Yc  data  points  to  plot  the  arbitrary 
cross  section  compression  surface.  We  recall  equations  (2- 
14)  and  (2-15)  modified. 

Vn«9  *  Ypo,  -  +crzsin$t  (2-42) 

Now  start  at  Y  and  march  away  in  small  increments  until  you 
reach  Y  .  At  each  step  along  the  way  stop  and  calculate  X 
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according  to  equation  (2-20)  and  R^  according  to  equation  (2- 

16).  We  can  now  determine  (R-J.,-:  using  (2-41).  With  ( R_.) 

cd  proj  3  CD  proj 

known  we  determine  X£,  Y£  with  the  following  simple 

tr igonometr ic  relations, 


Xc  -  X<RcbWR*b 


(2-43) 

(2-44) 


Yc  "  Y<RcbWR* 

Note  that  equations  (2-43)  and  (2-44)  are  similar  to 

equations  (2-19).  Also  in  equation  (2-41),  z^/1^  is  an 

arbitrary  input  value  depending  on  what  cross  section  is  being 

viewed.  For  example,  if  it  is  desired  to  view  the  cross- 

section  at  say  307.  from  the  leading  edge  then  z^/1^  ■  0.30  . 

We  now  have  a  systematic  way  to  view  the  waverider  shape 

at  any  arbitrary  cross  section  from  the  leading  edge  to  the 

base  plane.  Figure  2-9  below  gives  an  example  of  the  results 

using  this  procedure.  The  figure  shows  a  base  plane 

perspective  of  the  waverider  at  0.1  intervals  from  the  leading 

edge . 


Figure  2-9  Baseplane  Perspective  of  Arbitrary 

Crossec  tions 

2.5  The  Top  (planform)  View 

To  determine  the  planform  view  of  the  waverider  we 
consider  the  following  coordinate  system.  We  note  that  along 
the  freestream  surface  r&  is  a  constant,  and  if  we  project  to 
the  base  plane  we  see  that  r©  ■  r^  .  Therefore  we  have 

r©  *  rf(3  -  constant  (2-45) 

We  now  establish  a  relationship  between  1  and  1„  .  Me  can 
see  that  from  figure  2-10a, 

1  -  1„  -  2}!y«0)  (2-46) 

and  from  figure  2-10b,  9  *  R  when  r  is  drawn  to  the 

leading  edge  of  the  waverider  where  the  shock  intersects  the 
body.  Redrawing  figure  2-10b  we  have 
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Figure  2-10  Planform  and  Side  View  Coordinate  System 
( dimensional ) 


Now  using  basic  trigonometry  we  see  that, 

z,(y=0)  »  R016/tan|3  (2-47) 

but  using  small  angle  approximations  tanfl  ■  (3,  therefore 

zs(y-0)  »  Rfll  6/(3  *  R0l/a  (2-48) 

Substituting  (2-48)  into  (2-46)  we  have 

1  -  1.  “  Ro1/ff 

or  1^  =  1  (  l-RQl/cr)  (2-49) 


Figure  2-11  Side  View  (dimensional) 

Equation  (2-49)  gives  us  a  relationship  between  1  and 
Now  at  the  base  plane  r  *  1,  therefore 

1 ©^  -  rs(3  *  constant 

Dividing  through  by  8  and  rearranging  we  obtain 

9^/6  ■  rf|3/(18)  *  rfa/l 

Recall  that  ©  and  r  vary  for  a  given  shape  and  8,  13, 
do  not.  Now  9^/8  *  R^J#)  and  this  gives  us 

R#b($)  -  rf*/l 

For  small  angles  r}  *  zf  and  we  rewrite  (2-51)  as, 

z/1  *  Rfb(*)/a 

Now  recall  from  our  spherical  coordinate  system  that, 

y  **  r©sin$,  r  ■  r 
or  y  *  zGsinf 


(2-50) 
and  a 

(2-51) 

(2-52) 

(2-53' 

(2-54) 


Along  the  boundary  (where  the  shock  and  body  intersect)  © 
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3 .  Replacing  2  with  r$  and  dividing  through  by  1  we  have, 

y/1  =  rt3sin$/l  (2-55) 

but  r %/ 1  =  R^($)/a  from  (2-52).  Therefore  substituting 

yields,  1 

y/i  =  (2-56) 

Equations  (2-52)  and  (2-56)  give  us  a  way  of  drawing  the 
planform  view  normalized  with  1.  However,  it  may  be  more 
convenient  to  normalize  with  1^.  Using  (2-49)  to  accomplish 
this  we  have, 

V1,  *  -  R0>  (2-57) 

y,/  1  „  =  R^(#)6sin$/(l  -  R0/e)  (2-50) 

Equations  (2-57)  and  (2-58)  give  us  a  method  of  plotting 
our  waverider  as  viewed  from  the  top  normalized  with  1^.  The 
appendix  of  this  thesis  devotes  much  of  its  bulk  to  showing 
various  planform  and  baseplane  views  of  different  waveriders 
determined  by  a  choice  of  different  parameters  in  equation  (2- 
10)  . 

2.6  Three  Dimensional  Views 

Development  of  a  three-dimensional  view  of  an  arbitrary 
waverider  shape  can  be  an  extension  of  the  arbitrary  cross 
section  and  planform  views  presented  earlier.  Adamsli  has 
developed  a  code  using  this  analysis  to  present  a  3-D  view. 
Several  views  are  presented  in  appendix  H  using  this  code. 


CHAPTER  3 


Surface  Skin-Friction  Coefficients 

3.1  Introductory  Remarks 

In  this  chapter  we  wish  to  derive  expressions  for  the 
average  skin-friction  coefficients  for  both  the  freestream  and 
compression  surfaces.  Although  these  expressions  will  be 
approximate,  they  will  reflect  the  proper  trends  for 
variations  in  Reynolds  number,  Mach  number,  and  surface  wall 
temperature.  Since  the  pressure  is  constant  on  the  freestream 
surface  of  the  waverider  and  varies  slowly  along  an  inviscid 
streamline  on  the  compression  surface,  we  will  assume  that  the 
friction  coefficient  can  be  modeled  after  a  flat  plate  for 
the  freestream  surface  and  a  wedge  for  the  compression 
surface.  In  both  cases,  the  pressure  gradient  along  an 
inviscid  surface  streamline  is  zero.  We  shall  use  the 
momentum-integral  formulation  for  compressible  boundary 
1 ayers . 

3.2  Momentum- Integral  Equation  for  Flat  Plates  and  Wedges 
The  momentum-integral  equation  for  zero  pressure  gradient 

is  given  by  reference  18  as, 
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Where  9  is  the  compressible  momentum  thickness  given  by, 


(3-2) 


and  t  is  the  surface  stress  determined  by, 


3* 


(3-3) 


Here  x  and  y  are  measured  along  and  normal  to  the  surface  in 
the  usual  way,  as  shown  in  figure  3-1.  The  thickness  of  the 
boundary  layer  is  denoted  by  $*S(x).  The  subscripts  e  and  w 
denote  conditions  at  the  outer  edge  of  the  boundary  layer  and 
at  the  wall  surface. 


CondrtioflS 


Figure  3-1  Boundary  Layer  Definitions 
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To  account  for  compressibi 1 i ty  effects,  we  introduce  a 
transformed  normal  coordinate  Y  defined  by, 

dY  =  -Q-  dy  ( 3-4 ) 

He 

The  momentum  thickness  can  now  be  written  as, 


©■ 


(3-5) 


where  Y  *  A  is  the  location  of  the  outer  edge  of  the  boundary 
layer  for  the  transformed  coordinate.  The  wall  friction  is 
given  by, 


Tw 


(3-6) 


3.3  Approximate  Solution  Based  on  Assumed  Velocity  Profile 
An  approximate  solution  for  the  momentum  integral  equation 
can  be  obtained  by  using  an  assumed  velocity  profile  having 
the  form, 


(3-7) 


The  momentum  integral  and  wall  friction  can  now  be  expressed 
as 

0  =  A  a2 


( 3-Ba ) 
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where 


(3-Sb) 


(3-9a) 


(  3-9b  ) 


For  this  problem,  we  assume  that  the  pressure  gradient  is 
zero  and  that  the  wall  temperature  is  a  constant,  and  thus 
that  ctj  and  a2  are  constants  that  depend  only  on  the  assumed 
velocity  profile.  The  momentum-integral  equation  (3-1)  can 
now  be  expressed  as, 


AdA  _  _  «1 

dx  ~  w^eue  <*2 


(3-10) 


where 


Cw 


Pm  H. 


Pa  M. 

0  6 


(3-11) 


For  this  problem,  the  right  side  of  equation  (3-10)  is  a 
constant.  Thus  we  can  integrate  it  without  difficulty, 
subject  to  the  boundary  condition  A(0)  =  0,  and  obtain 


A 


< 

N 


ia 

on 


ifaA. 


(  3-12) 
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We  thus  obtain  the  well-known  result  for  flat  plates  that  the 
transformed  boundary- 1 ayer  thickness  grows  like  the  square 
root  of  x  . 

The  local  skin  friction  can  now  be  obtained  from  equation 
(3-Qb).  The  local  sk in-f r ic tion  coefficient,  based  on  the 
external-flow  conditions  is, 


where 


yj  2  «1  «2~Cw 


s  PPUft  X 


(3-13) 


(3-14) 


is  the  Reynolds  number  based  on  the  external  flow  conditions. 
Equation  (3-13)  has  the  same  form  as  that  obtained  from  the 
classical  Blasius  equation  for  constant-density  flow  past  a 
fiat  plate.  For  compressible  flow,  the  factor  CB  depends  on 
the  assumed  viscosi ty- temperature  relation.  The  factor 
depends  on  the  assumed  velocity  profile. 

3.4  Viscosi ty-Temperature  Power  Law  Relation 

In  this  analysis,  we  shall  assume  the  viscosity- 
temperature  power-law  relation 


(3-15) 
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where  u>  is  a  constant  usually  taken  to  lie  in  the  range 
.5  1  u>  <  1 ,  depending  on  the  temperatures  involved  in  the 

problem.  Generally,  U  is  taken  to  be  about  3/4. 

For  the  special  case  OJ  =  1 ,  equation  (3-15)  takes  the 
form 


Js.  JL  =  i 

T  Pe 


(3-16) 


Since  the  pressure  is  a  constant  across  the  boundary  layer, 

we  have  for  a  thermally  perfect  gas  (P  =  PRT)  that  -I?-  =  ^  , 

T 

and  thus  that 


P  P 


(3-17) 


For  this  case,  therefore,  we  have  that  C  «  1. 
1,  C|f  from  equation  (3-13)  has  the  same 
incompressible  flow,  provided  that  f>  and  p 
at  the  external  flow  conditions. 

If  u  *  1  ,  equation  (3-15)  can  be  written 


Thus  for  u>  = 
form  as  for 
are  evaluated 

in  the  form 


where 


c 


JLtt 

p  p 

e  e 


(  3-18) 


(3-19) 
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Equation  (3-16)  is  sometimes  referred  to  as  the  Chapman- 
Rubesin  form  of  the  viscosi  ty- temperature  law,  when  C  is  taken 
as  some  constant  or  some  function  of  x,  but  not  as  a  function 
of  y.  For  our  purposes  here,  since  T  is  a  function  of  y,  we 
will  consider  that  C  is  also  a  function  of  y,  given  by 
equation  (3-19).  It  corresponding  1 y  follows  that, 


Cw 


(3-20) 


3.5  Assumed  Velocity  Profile 

We  shall  assume  that  U/Ue  =  f(Y/A  )  has  the  specific  form 
following  the  Karman-Pohl hausen  technique.  Namely,  a  fourth- 
order  polynomial  having  the  form, 


U 

Ue 


(1-A-B-C) 


(3-21) 


This  expression  satisfies  the  condition  at  the  outer  edge  of 
the  boundary  layer  that  U/Uf(l)  *  1.  The  constants  A,  B,  and 
C  are  as  follows  to  meet  boundary  conditions. 


A 


9  ( U/Ue )  1  .  ai 
a(Y/A  )JW 


B  = 


2 


3  (U/U«) 

.oUTlT)2 


(  3-22a ) 


(  3-22b  ) 
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C  =  (-2  -  B)  ( 3-22c ) 

To  specify  A  and  B  and  thus  C,  we  start  by  noting  from  the 
boundary- 1 ayer  momentum  differential  equation  that  at  the  wall 


(3-23a) 


provided 
terms  of 


the  pressure  gradient  along  the  wall  is  zero, 
the  transformed  variable  V,  this  is 


In 


( 3-23b ) 


Since  Me  does  not  depend  on  Y,  we  can  expand  equation  (3-23b) 


and  get 


(3-24) 


Equation  (3-22)  provides  a  relation  between  the  first  and 
second  derivative  of  the  velocity  at  the  wall,  that,  is, 
between  A  and  B.  It  is  useful  to  place  it  in  an  alternative 
form.  We  note  along  with  the  expressions  (3-3)  and  (3-6)  at 
the  wall  that  we  have  similar  expressions  for  the  heat  flux 


at  the  wall 


q 


w 


(3-25) 


Consequently,  equation  (3-24)  can  be  written  alternatively  as 


/  a2u> 

k  (u- 1) 

JL!». 

/  d  u  \2 

'»■  T* 

Cp  Tw 

(3-26) 


where 

pr  “  fyS/k.  (3-27) 

is  the  Prandtl  number  evaluated  at  the  wall,  and  Cp  is  the 
specific  heat  at  constant  pressure.  Me  assume  that  both  Pf 
and  Cp  are  constants. 

It  is  further  useful  to  introduce  the  Stanton  number  as 
the  dimensionless  heat-transf er  coefficient* 

3  _  -  (3-28) 

^  U#  Cp  (  Tw  -  Tflw) 

Here  the  Stanton  number  is  based  on  the  external  flow 
conditions,  and  TM  is  the  adiabatic  wall  temperature,  to  be 
specified  later.  Note  that  if  heat  is  transferred  from  the 
fluid  to  the  wall,  then  and  qw  <  0.  In  terms  of  Sj  and 

Cf,  defined  by  equation  (3-13),  we  can  rewrite  equation  (3-24) 
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or 

B  -  (t*>~  Tow  )  2  Pf  ST  >2  ( 3-29b ) 
2  Tw  Cf 

The  factor  2ST/Cf  is  known  as  the  Reynolds-analogy  factor. 
It  will  be  specified  later.  Note  when  u  *  1  that  B  »  0,  and 
also  when  C  is  a  constant  or  a  function  of  x  only  that  B  again 
is  zero. 

Since  8  is  a  function  of  A,  we  need  a  condition  on  the 
assumed  profile  (3-19)  to  determine  A.  To  do  this,  we  assume 
the  following  condition  at  the  outer  edge  of  the  boundary 
1  ayer : 


ra(U/U«)  I 

Jr.. 


This  gives, 

A  +  2B  3C  +  4 ( 1 -A-B-C )  -  0 
substituting  the  relations  for  B  and  C  yield, 

DA2  +  3A  -  6  -  0  ( 3-31 ) 


where 

D  e  ( u  Z  ^  Tqw  ^  2 

Tw  Cf 


2 


(3-32) 
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Solving  the  quadratic  equation  for  A  gives 

A  =  12/(3  ♦  T(9+24D))  (3-33) 

When  D  *  0,  then  A  *  2,  and  the  velocity  profile  reduces  to 

2  ft) 

3.6  Effects  of  Mach  Number  and  Heat  Transfer  on  the  Friction 
Coefficient 

The  momentum- thic kness  parameter  a7  can  now  be  evaluated 
from  the  basic  definition  equation  (3-9a)j 

a7  -  1/35  +4A/15  -  A2/9  +  DA2(1/14  -9A/140  -  DA2/105)  (3-35a) 

Since  A  depends  on  the  factor  D  from  equation  (3-33),  then  c«2 
depends  on  D  only.  Using  the  relationships  from  reference  17, 

A  ■  2  +  •  a  -  -  2  B  (3-35b) 

o 

and  that  D  *  BA  allows  us  to  write  (3-35a)  in  the  same  manner 
as  is  done  in  the  same  reference  namely, 

a2-  JL  -  ±\ 

63  \  5  15  144/ 


(  3-35c  ) 
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When  D  =  0,  that  is,  when  w  =  1  or  when  (no  heat 
transfer),  then  A  =  2 .  We  then  have  »  0.1175.  Likewise 
for  this  case  we  have  C#  ■  1 ,  and  the  factor  appearing 
in  the  friction  coefficient  C{e,  equation  (3-13),  takes  the 
va  1  ue , 

/  (  2 o jQ j ~  0 . 685  . 

The  correct  value  for  this  factor  evaluated  by  the  Blasius 
profile  form  laminar  boundary- 1 ayer  theory  is  0.664.  Thus, 
for  this  case,  the  profile  from  equation  (3-34)  leads  to  an 
error  of  3  percent. 

Although  this  error  is  quite  acceptable  for  our  present 
purposes,  the  important  utility  of  the  present  analysis  lies 
in  the  result  for  equation  (3-35),  that  produces  combined 
results  for  u)  #  i  and  for  wall  heat  transfer,  and  for  which 
there  are  no  corresponding  analytical  results  from  the 
di f f erential-equation  analysis  of  compressible  laminar 
boundary- 1 ayer  theory.  To  proceed  further,  we  need  to 
evaluate  the  factor  D,  and  to  do  this  we  need  expressions  for 
the  adiabatic  wall  temperature  and  the  Reynolds  analogy 
factor.  We  shall  thus  assume  the  relations 


'ow 


1+  v|Pr 


i~P~  7-  \  M 


e 


2ST/C, 


P 


•2/3 


(  3-36) 


(3-37) 
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which  are  well  known  approximations  for  flat  plates  and 
wedges.  The  factor  D  thus  becomes 


(3-38) 


With  these  results,  we  now  have  the  general  functional 


dependence 


Cfe"  f  ^e***^.  t*  >pr.  y) 


(3-39) 


By  means  of  equation  (3-20)  we  can  write  equation  (3-13)  as 


(3-40) 


with  the  factor  4"  ( 2cij  aj)  being  a  function  of  D  as  given  by 
equation  (3-38).  For  adiabatic  walls,  the  results  simplify 
since  D  ■  0,  and  we  have 


The  factor  Cf#T(RM#)  is  plotted  versus  Mf  in  figure  3-2  for 
1  *  1*4,  w  »  0.7,  Pf  ■  0.72.  It  can  be  seen  that 


increasing  the  ex terna 1 -f 1 ow  Mach  number  can  reduce  the  local 
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skin-friction  coefficient  considerably. 


Figure  3-2 

Figure  3-2  also  shows  two  curves  for  when  the  wall 

temperature  is  held  fixed,  that  is,  for  T  /Tf  ■  1  and  6.  Again 

these  curves  are  calculated  by  means  of  equations  (3-39),  (3- 

36),  (3-33)  and  (3-31).  Again,  the  local  skin  friction 

coefficient  decreases  as  M.  increases.  These  curves  are  in 

c 

agreement  with  corresponding  results  obtained  by  numerical 
integration  of  the  similarity  boundary-layer  equations  of 
reference  19  and  20. 


3.7  Application  to  Waverider  Freestream  Surface 
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For  -flat  plates  the  external  conditions  are  the  same  as 
the  freestream  conditions.  Since  the  inviscid  surface 
streamlines  are  straight  lines  with  zero  pressure  gradient  on 
the  upper  freestream  surface  of  the  waverider,  we  can  apply 
the  previous  formulas  for  Cfe  by  replacing  the  subscript  e  with 
the  subscript  ®.  In  doing  this  we  ignore  the  fact  that  there 
is  actually  a  cross  flow  in  the  boundary  layer.  This  occurs 
because  the  boundary  layer  originates  at  different  locations 
along  the  leading  edge.  In  a  subsequent  section,  we  will 
obtain  an  average  skin  friction  coefficient  for  the  freestream 
wetted  surface  by  integrating  the  local  skin-friction 
coefficient  along  each  straight  streamline  on  the  freestream 
surface  and  then  integrating  again  across  the  span  of  the 
surface.  This  approximate  result  can  be  obtained  straight 
forwardly  and  provides  a  simple  means  for  making  parametric 
studies  of  friction  drag. 

3.8  Compression  Surface  of  the  Waverider 

For  the  compression  surface  of  the  waverider,  the  external 
conditions  for  the  boundary  layer  are  not  the  same  as  the 
freestream  conditions  because  the  inviscid  surface  streamlines 
pass  through  the  conical  bow  shock.  To  express  the  local 
skin-friction  coefficient  and  the  local  Reynolds  number  in 
terms  of  the  freestream  conditions,  we  therefore  rewrite 
equation  (3-13)  as  follows: 


(3-42) 


where 


P„  U»  * 

R«*- =  H. 


For  a  thermally  perfect  gas  we  have 


(3-43) 


ik  a  Jl_  (3-44) 

^cc  Peo  Te 


Thus,  when  the  density  is  eliminated  in  favor  of  the  pressure, 
equation  (3-42)  becomes 


2 


(  3-45) 


The  inviscid  surface  streamlines  on  the  compression  side 
of  the  waverider  lie  in  a  plane,  and  the  pressure  increases 
slowly  along  the  streamline.  Although  the  boundary  layer  has 
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a  cross  flow  and  is  thus  actually  three  dimensional,  we  assume 
these  effects  are  small  enough  to  ignore,  and  thus  the 
boundary  layer  along  an  inviscid  surface  streamline  behaves 
nearly  as  that  of  a  flat  plate.  further,  we  assume  that  the 
external  flow  conditions  are  given  by  the  conditions 
immediately  behind  the  conical  bow  shock.  If  we  evaluate 
these  conditions  by  hypersonic  small-disturbance  theory  then 
we  have 


Ue 


U 


*  1 


(  3-46a) 


JL-  «  1  +  7  Kg  (3-46b) 

P  ° 

'ao 


+  2  =  Te_ 
+  2  T» 


1  + 


(7-1) 

2 


(  3-46C  ) 


where 


(7+1) 


(  3-46d  ) 


and  Kg  =  .  Equation  (3-45)  now  takes  the  form 


yj  2  a, Q2 


(3-47) 


The  factor  ( ot  jOtj  ^  e  must  be  evaluated  such  that  and  Me  in  the 
factor  D  are  calculated  in  terms  of  the  infinity  conditions 
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by  means  of  equation  (3-46c). 

Ule  may  simplify  equation  (3-47)  by  writing  it  in  terms  of 
the  freestream  skin  friction  coefficient  given  by  equation  (3- 
40).  Doing  so  gives, 


K(C4  ) 


foo'fs 


(3-48) 


where 


(3-49) 


In  the  factor  K  the  dominant  term  is  ■T(Pe/'P#)  since 
S  (  ( a^a^)  ( a^)  )  is  close  to  unity.  In  fact  it  is  unity  when 
Pr  -  1. 

3.9  Average  Skin  Friction  Coefficients 

The  average  skin  friction  coefficient  is  determined  by 
summing  up  all  of  the  local  coefficients  along  each  streamline 
and  then  dividing  through  by  the  wetted  area  of  the  surface. 

Considering  figure  3-3  below,  we  have  for  the  average  skin 
friction  coefficient  for  the  freestream  surface 


C 


f  f 


(3-50) 


dz  ds 
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In  this  section  we  will  skip  the  details  of  the  analysis  of 
the  above  equation  and  save  it  for  appendix  A.  However,  after 
this  analysis  we  have, 


-  "«/<7  Cf#  2^1- VWl  +  CdX/dY?  dY 


'/l  +  (<JX/dyj,<1Y 


(  3-51  ) 


The  reader  may  choose  any  numerical  integration  scheme  desired 
to  evaluate  the  integrals. 

Now,  to  determine  the  average  skin  friction  coefficient 
for  the  compression  surface  we  will  make  the  assumption  that 


the  wetted  area  for  the  compression  surface  is  approximately 
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the  same  as  that  of  the  freestream  surface.  Using  this 
approximation  gives  for  the  compression  surface, 


c  2  pi-  Rt/»  'J 1  +  (dX/<fY  j  dY 
JO“Rt/«  )  +  (dX/dY  j  dY 


(3-52) 


In  equations  (3-51)  and  (3-52)  above  and  C<c  are  the 
local  skin  friction  coefficients  determined  using  the  analysis 
of  the  previous  sections. 

3.10  Estimation  of  Total  Drag 

A  good  discussion  on  various  types  of  Drag  is  presented  in 
reference  21.  Here,  we  will  begin  with  the  following 
definition 


D  *  +  Db  (3-53) 

where  is  the  wave  drag,  Dj  is  the  skin  friction  drag  and 
is  the  base  drag. 

The  wave  drag  exists  at  supersonic/hypersonic  speeds  and 
is  caused  by  static  pressure  gradients  across  the  shock  wave. 
We  will  define  the  wave  drag  as 

D»  =  qAbC#  (3-54) 

2 

where  q  =  1/2  p#U  A^  is  the  base  area  of  the  waverider  and 
is  the  wave  drag  coefficient. 
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The  skin  friction  drag  is  caused  by  shearing  stresses  over 
the  wetted  surface  of  the  body.  We  will  define  the  friction 
drag  as 

Df  =  9^ft(Sirf  +  S«c)  (3-55) 

where  and  represent  the  wetted  surface  area  of  the 

freestream  and  compression  surfaces  respectivel  y .  Also,  Cft 
is  the  average  total  skin  friction  coefficient.  However, 
using  the  analysis  of  the  previous  sections  we  can  rewrite  ( 3- 
55 )  as 

Df  =  q(S-fC7f  +  S^C^)  (3-56) 

or  using  the  results  of  equation  (3-48)  we  have, 

D,  *  OC^S.,  .  KS„>  (3-57) 

The  base  drag  is  caused  by  the  blunt  end  of  a  body.  This 
is  because  the  pressure  at  the  base  is  usually  not  the  same 
as  the  surrounding  infinity  conditions.  We  will  ignore  this 
effect  however,  as  is  usually  done  when  comparing  different 
configurations^.  Thereby  setting  *  P#. 

Ignoring  the  base  drag  and  combining  equations  (3-51)  and 
(3-50)  into  (3-49)  gives 

D  =  q  ( C^Ajj  +  +  C<cSc) 

or  using  the  results  of  (3-57)  we  have 

D  =  q[C(|Ab  +  +  KS^)] 

It  may  be  more  convenient  to  have  the  total  drag  in  terms 
of  the  average  total  skin  friction  coefficient.  Equating  (3- 


(3-58) 


(3-59) 


55)  and  (3-57)  and  solving  for  yields, 

-  Cfit  <  S,1  *  V/(S.(*  KS.)]  (3-60) 

Substituting  (3-60)  into  (3-59)  gives, 

o  -  qCCjft,  -  c7t(Srt  +  S.n  (3-61) 

where 

cTt  -  ^t‘S„  *  KS,)/(Slri  -6  S„)J  (3-62) 

In  the  appendix  to  this  thesis  we  present  the  lift-to-drag 
ratio  as  well  as  the  average  total  skin  -friction  coefficient 
for  each  particular  configuration  shown. 


CHAPTER  4 


Comparing  Configurations 

4.1  Introductory  Remarks 

Here  we  will  describe  what  happens  as  we  change  different 
parameters  of  the  waverider  configuration.  In  order  to  do 
this  we  must  establish  a  basic  shape  from  which  we  will  base 
our  comparisons.  Once  the  basic  configuration  is  chosen  then 
we  will  discuss  what  happens  when  we  change  the  coefficients 
of  Y  of  equation  (2-10)  namely,  bj,  and  bfc.  After  this  we 

will  compare  shapes  with  different  values  for  Rgl  and  6  and 
include  a  discussion  on  the  Mach  number  and  Reynolds  number. 

4.2  The  Basic  Waverider  Conf iguration 

The  basic  shape  is  chosen  to  have  the  parameters  shown  below 
in  table  4.1 

Table  4-1 

Basic  Configuration  Parameters 


* 

<5 

Mach 

Number 

Reynolds 

Number 

7 

50° 

8 

10 

1X106 

1.4 

*> 

b2 

b4 

slope 

0.5 

0.903 

-1.0 

0.437 

35° 

46 
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The  basic  parameters  were  chosen  arbitrarily  except  for  bj 
and  b4>  The  value  of  b^  was  chosen  such  that  the  slope  at  the 
edge  of  the  waverider  freestream  surface  was  25*.  This  value 
is  in  between  the  limits  imposed  by  condition  iii)  of  chapter 
2.  The  value  of  bj  was  fixed  according  to  equation  (2-13). 
These  basic  configuration  parameters  resulted  in  the  waverider 
conf iguration  shown  below  in  figure  4-1. 


Figure  4-1  Basic  Configuration  Base  Plane  and  Planform 
Views 

We  will  now  vary  these  parameters  and  see  how  the  basic 
conf iguration  changes. 


2.3  Changing  the  Coefficients  of  Y  of  the  Equation  Describing 
the  Freestream  Edge  in  the  Base  Plane  With  Consideration  of 
the  Edge  Slope 
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Here,  we  will  hold  the  parameters ,  #j,  6,  M ,  Rey,  and 

constant  and  vary  in  a  systematic  way  the  coefficients  of  Y 
equation  (2-10). 

Suppose  we  want  to  change  bj  from  a  -1.0  as  shown  in  table 
4.1,  to  a  -2.0.  This  is  easily  done.  However,  due  to  equation 
(2-13)  changing  b4  will  result  in  a  different  value  of  bj  as 
well.  We  must  also  keep  in  mind  the  imposed  condition  of  the 
slope  at  the  edge  of  the  waverider.  Changing  b^  will  also 
change  the  slope  at  the  edge  of  the  waverider.  We  can 
however,  keep  the  slope  at  a  constant  value  by  letting  bfc 
change  by  a  suitable  amount.  Table  4-2  below  shows  what 
happens  when  we  change  bj  from  -1.0  to  -2.0,  holding  RQ  and  the 
edge  slope  constant.  Also,  figure  4-2  shows  the 

corresponding  waverider  configuration. 

Table  4-2 


$ 

<5 

Mach 

Number 

Reynolds 

Number 

y 

50° 

8 

10 

1X106 

1,4 

Ro 

b2 

k>4 

*>6 

slope 

0.5 

1.405 

-2.0 

0.935 

. 

25° 

Figure  4-2  Effect  of  Changing  b4 


This  discussion  shows  the  inherent  complexities  of  changing 
what  may  seem  to  be  a  rather  simple  variable.  By  changing  b^ 
we  discovered  that  bj  and  b&  must  also  change. 

Suppose  now  instead  of  changing  bj  we  want  to  change  the 
freestream  edge  slope  from  25*  to  its  maximum  limit  of  50*. 
This  is  accomplished  by  changing  the  value  of  b^  or  b^. 
However,  for  the  sake  of  simplicity  we  will  hold  b^  at  its 
previous  value  of  -1.0.  Setting  dX/dY  *  50  and  solving  gives 
a  new  value  of  b^  and  bj.  Table  4-3  below  shows  this,  and 
figure  4-3  shows  the  correspond ing  waverider  configuration. 

A  similar  approach  to  the  above  analysis  can  be  done  if  we 
want  the  slope  of  the  freestream  edge  at  its  minimum  value. 

We  must  also  ask  ourselves  what  are  the  aerodynamic  and 
geometric  consequences  of  changing  the  values  of  b^.  Analysis 
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Table  4-3 


* 

5 

Mach 

Number 

Reynolds 

Number 

7 

50° 

8 

10 

1X10* 

1.4 

Ro 

b< 

mm 

slope 

0,5 

0.795 

1.0 

0.544 

ft 

50 

has  shown  that  the  viscous  lift  to  drag  ratio  may  increase 
or  decrease  with  a  change  in  b^.  The  deciding  factor  is  the 
initial  value  of  b^.  However,  the  change  in  magnitude  of  the 
viscous  L/D  appears  to  change  only  rather  slightly. 


Figure  4-3  Effect  of  Changing  Freestream  Edge  Slope 


Similarly,  the  effects  of  changing  the  edge  slope  in  general 
has  a  minor  effect  on  the  aerodynamic  values.  However,  a 
slope  of  SO*  usually  gave  higher  L/D  ratios;  although  in  some 
cases  this  was  not  true. 

4.4  Changing  Rfl  in  the  Equation  Describing  the  Freestream 
Edge 

When  changing  Rfl,  we  might  say  that  there  is  no  effect  on 
the  edge  slope.  However,  a  further  investigation  tells  us 
that,  even  though  dX/dY  has  no  RQ  terms  in  it,  the  value  of  bj 
does  depend  on  RQ.  Since  the  slope  is  dependent  on  b2  we  see 
that  changing  Ro  will  result  in  a  changing  tip  slope.  We  can 
however  keep  the  tip  slope  constant  by  changing  b^  as  done 
before . 

A  more  dramatic  affect  caused  by  changing  Rfl  is  the  non- 

dimensional  length  of  the  waverider.  This  is  easily  seen  by 

equation  (2-38).  Also,  figure  (4-4)  below  shows  this  result. 

We  can  see  that  small  values  of  R„  will  result  in  more  slender 

o 

configurations.  We  define  the  aspect  ratio  for  these 
con f igurations  as  b/1^  Where  b  is  the  maximum  distance 
across  the  base  of  the  waverider.  Therefore  the  smaller  the 
aspect  ratio  the  more  slender  the  waverider.  Aspect  ratios 
are  given  for  each  configuration  shown  in  the  Appendix. 


Figure  4-4  Effects  of  Changing  Rfl 

Analysis  shows  in  general  that  small  values  of  RQ  usually 

result  in  better  viscous  lift  to  drag  ratios. 

Another  effect  on  the  waverider  that  occurs  with  a  change 

in  is  the  shock  stand  off  distance.  Smaller  values  of  R_ 
o  o 

result  in  the  shock  being  a  greater  distance  away  from  the 
compression  surface.  This  can  be  seen  in  the  above  figure. 
There  are  other  factors  that  determine  the  shock  stand  off 
distance  and  these  will  be  discussed  as  they  occur. 

Changing  RQ  also  significantly  changes  the  shape  of  the  base 
plane  configuration  as  is  shown  below  in  figures  4-5a,b. 

In  figures  4-5a,b  the  value  of  the  edge  slope  and  the  value 
of  b^  is  kept  constant.  In  general,  larger  values  of  RQ  cause 
the  freestream  surface  to  be  flatter. 


Figure  4-5a  RQ=0.7,  b^-l.O,  Figure  4-5b  Basic  Config 

slope=25°  Rq=0.5,  b^-1.0,  slope=25* 


4.5  Changing  the  Parameter 

When  we  change  <fcj  the  most  significant  change  in  the 
waverider  shape  that  is  first  noticed  is  the  change  in  maximum 
span  or  width;  all  other  parameters  remaining  the  same. 
Equations  (2-14)  and  (2-15)  give  us  the  positive  and  negative 
most  values  of  Y,  which  allow  us  to  determine  the  span  of  the 
waverider.  Clearly  a  decrease  in  decreases  the  span.  A 
decrease  in  the  span  decreases  the  aspect  ratio  and  increases 
the  slenderness  of  the  waverider. 

If  we  change  holding  and  b^  constant,  we  get  a  new 
value  for  b2  according  to  equation  (2-13).  Also,  b6  will  have 
to  change  to  keep  the  slope  at  the  edge  constant  or  within 
limits  according  to  condition  iii).  Figures  4-6a , b  show  two 
configurations  with  identical  8,  b^,  and  tip  slope  but  with 
different  '  s 


Me  see  from  this  that  a  change  in  causes  a  significant 
change  in  the  waverider  shape;  not  only  in  span  but  also  in 
the  curvature  of  the  freestream  and  compression  surface. 
However,  there  is  no  change  in  the  length  of  the  waverider. 


4.6  Changing  the  Parameter  S 
When  we  change  6  the  value  of 
on  the  standoff  distance  between 
surface  of  the  waverider  body, 
definition  of  a  from  chapter  3. 

a2  .  12112 

2 

where  Kj  =  M#S . 

We  can  see  from  equation  (4-. 


changes.  This  has  an  effect 
the  shock  and  the  compression 
To  see  this  we  look  at  the 


)  that  a  decrease  in  6  will 


result  in  an  increase  in  a. 


This  causes  the  separation 
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between  (3  and  6  to  increase;  which  affects  the  shock  standoff  j 

distance.  The  standoff  distance  becomes  important  when  ] 

considering  design  of  engine  and  engine  inlet  conditions.  j 

1 

Below  in  figures  4-7a,b  we  show  different  waveriders  with  the  3 

same  Rfl,  $j,  b4 ,  and  tip  slope  but  with  different  6's. 


i 

l 

I 

i 

i 

) 


Figure  4-7a  R0=O.5,  b^=-l .  0  , 

6=6,  $j=50,  slope=25* 

Again,  changing  6  and  therefore  <7,  results  in  a  change  of 
bj  according  to  equation  (2-13)  and  a  change  in  b4  or  b&  to 
meet  the  slope  condition  on  the  edge  of  the  waverider 
freestream  surface.  These  changes  cause  the  curvature  of  the 
freestream  and  compression  surfaces  to  vary. 

4.7  Changing  the  Mach  Number 

Varying  the  mach  number  has  a  similar  effect  as  changing  6 
as  we  see  from  equation  (4-1).  An  increase  in  M  will  cause 


R=0.5,  b^=-l  .  0 ,  S»8,  4>j“50, 

si ope=25 • 
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the  standoff  distance  to  be  smaller  just  as  an  increase  in  6 
will  do.  The  only  difference  will  be  the  magnitude  of  the 
change.  Again  we  must  take  into  account  how  changes  affect 
cr  and  therefore  bj*  b^,  and  b^. 

4.0  Changing  the  Reynolds  Number 

We  realize  that  M  has  a  role  in  the  determination  of  the 
Reynolds  number.  However,  lets  assume  the  Mach  number  will 
remain  constant.  This  leaves  alt i  ude  and  the  length  of  the 
waverider  as  the  deciding  factors  c  the  value  of  the  Reynolds 
number . 

Figure  4-8  below  shows  L/D  vs.  Reynolds  number  for  the  basic 
configuration  shown  in  figure  4-1. 


L/D  va  P  Holding  Mach  constant 


Figure  4-8  L/D  vs.  R  holding  Mach  constant 


CHAPTER  5 


Conclusion 


5.1  Summary 

The  two  primary  goa’  established  in  chapter  1  have  been 
accomplished.  In  the  appendices,  a  catalogue  of  over  100 
shapes  is  presented.  This  catalogue  shows  variations  of 
waveriders  that  can  be  generated  by  axiaymmetric  flow  past  a 
cone  and  is  the  result  of  tne  first  goal. 

As  discussed  earlier  there  can  be  an  infinite  number  of 
waverider  shapes  so  as  yet  the  catalogue  is  incomplete. 
However,  one  may  get  an  idea  and  a  feel  for  waverider  shapes 
by  studying  the  appendices. 

Chapter  3  is  the  result  of  the  second  primary  goal.  Here 
we  developed  an  analytical  method  of  estimating  the  skin 
friction  coefficient  --  where  as  before  only  numerical  methods 
existed.  An  average  total  skin  friction  coefficient  is 
presented  in  the  appendices  for  each  arbitrary  waverider 
shape . 


5.2  Resul ts  and  Discussions 

The  best  waverider  shape  depends  on  the  designers  main 
goals.  As  in  all  design  processes  one  is  usually  faced  with 
compromise.  For  example  the  waverider  with  the  best  lift-to- 
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drag  ratio  will  probably  not  have  the  best  volume-to-area 
ratio . 

The  following  table  is  a  tool  to  help  in  the  design 
process.  In  the  table  we  can  see  the  effects  of  certain 
geometric  and  aerodynamic  values  when  we  vary  parameters  that 
affect  the  shape  of  the  waverider.  When  using  the  table  the 
results  are  given  such  that  all  other  parameters  are  held 
constant.  For  example,  lc  ing  will  in  general  result  in 
a  lower  value  of  L/D  the  other  parameters  are  held 

constant;  namely,  RQ  and  b^.  The  table  doesn't  give 
magnitudes  of  change  only  direction  of  change.  In  some  cases 
the  effect  of  a  change  cannot  be  determined  ahead  of  time. 

An  example  would  be  the  effect  in  the  volume-to-area  ratio 
with  a  change  in  b^.  Lowering  b4  may  result  in  either  an 
increase  or  decrease  in  the  volume-to-area  ratio  depending  or. 
the  original  value  of  b^. 

TABLE  5-1 

EFFECT  OF  LOWERING  (*,,  RQ,  and  b4) 
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Following  this  chapter  in  appendices  B  through  H  a 
catalogue  of  over  100  shapes  is  presented.  All  of  the 
configurations  presented  are  based  on  Mach  number  of  10  and 
Reynolds  number  of  10^.  The  lift-to-drag  ratios  were 
calculated  and  found  to  lie  in  a  range  from  4.32  to  5.06  while 
the  average  skin  friction  coefficient  ranged  from  0.00176  to 
0.00205.  The  volume  to  area  ratio  V^/S^  ranged  from  0.224  to 
0.153.  In  the  table  under  each  view  are  various  geometric  and 
aerodynamic  values.  Some  of  them  are  compared  with  what  we 
call  the  idealized  waverider.  Reference  14  describes  this 
ideal  waverider.  The  freestream  surface  of  the  ideal 
waverider  is  given  by, 

X  *  |Vcot*j|  (  5-1  ) 
where  the  lines  about  the  term  on  the  right  represent  the 
absolute  value.  Note  than  in  the  above  equation  there  is  no 
term  representing  the  value  of  RQ.  This  is  because  Rfl  is  zero 


in  this  case. 
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Appendix  A 

Development  of  Equation  3-51 


We  begin  by  establishing  a  few  definitions  that  will  be 
used  later  in  this  analysis.  Consider  the  second  part  of 
figure  3-3  and  we  may  say  without  much  effort, 


ds  "  ^(dx)V(dy)* 

“  >/l  +  (dx/dy  f  dy 


Next,  consider  equations  (2-5)  and  (2-20).  From  these  we  can 
show  that, 

dx/dy  *  dX/dY  A. 2 

and  also  that, 

dx  *  1 SdX ,  dy  =  1 SdY  A. 3 

Now  considering  the  first  part  of  figure  3-3  we  may  write 
that  the  wetted  surface  area  of  the  freestream  surface  is, 


S 


W  to  ” 


A. 4 


Now  since  we  are  using  small  angles  we  have, 


Z5  "  r! 


There  f ore , 


A.  5 


A .  6 
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to  give, 


Now  use  equation  (2-29)  and  that  * 


Wee 


2  e  J(  i  -  V<0ds 


A. 7 


and  substituting  A.l  yields, 


ry‘  , _ T 

SWco“  22  l(  1  -  Rt/tr)  J  1  +  (d  x/dyf  d  y 

*  e\ 


A.0 


Now  if  we  substitue  relationships  A. 2  and  A. 3  we  have, 


I 

S  -  2  t*a  |Cl~  >/ 1  (dX/dY  f  dY 

Woo  J  a 


A. 9 


Recall  equation  (3-50) 


Cf,  .  -t-* 


J  J  cf-d<  d8 


JJ  dzd8 


A.  10 


where  C<#  is  the  local  skin  friction  coefficient  for  the 
fresstream  surface.  Substituting  A. 4  into  A, 10  yields, 


f  f 


sf  J/' 

Woo  i  it 


C.  dz  ds 

T  oo 


A.  11 


Using  equation  (3-40)  and  replacing  the  subscript  e  with  » 
according  to  section  3.7  yields. 
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,f  sw.  Jxi  V  Rexe  'Tw' 


ds 


=  _i_  v[2~^  llzY*~  f  f*  di— 

Sw«  ^Tw  *  >.n/  2  -  Z, 


>/2 


a,06 


s— 

4  \J  2 

SWoo  JPJLMl  XTW 

JL_  >d 


ds 


«S 


(r )  2  '7  ^ 


z«a  ds 


*/f 


S"~  'Jfi0'l,ew/Mb0 

Now  recall  from  equation  (3-28)  that, 


ds 


1„  =  1(1  -  Rq/<7) 


substituting  yields, 


—  .  JL  !lH5l  J ^ 

"  s-  WJ  J.  r‘ 


ds 


A. 12 


A.  13 


A.  14 


We  may  now  make  the  substitutions  of  equation  (3-40)  and 
equations  A.l,  A. 2,  and  A. 3  to  yield 


_  ^  ^"Vo  f>/i-Rt/<T  ^  1  +  (dX/dY  f  dY 


a.  i; 


APPENDIX  B 


Presented  here  are  configurations  that  have  the  following 

parameters  held  fixed: 

=  60,  6*8,  M  *  10,  Rey  *  1X106 

The  following  parameters  are  varied: 

R„  from  0.2  to  0.4 
o 

from  0.0  to  -3.0 

The  values  of  bj  and  b&  are  such  that  equation  2-13  and 
condition  (iii)  holds. 
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Figure  B.l 


Phi  -60.0 


CftaT  : 
S»/Sp  : 
ib/ibi  s 
b/li  = 


1 

'able  B.l 

m 

Hacb  =10  0 

8ej  s  0.101+07 

6aua  =  1.40 

I  :  0.200  ♦  1.057  Y*2  ♦  -1.000  T‘4  ♦  0.352  T‘6 


0.001947 

2.152777 

0.979190 

0.372399 


Cl  s 


7"(2/3)/Sp  : 
¥/?l  = 


SSD  s 


0.043795 

0.212266 

0.774453 

0.300503 


CO  2 

( L/D ) t Is  2 
Sp/Spi  = 


0.009691 

4.519036 

0.901519 


Figure  B. 2 


Figure  B.3 


Table  B.3 


Phi  -60.0  Delta  =  8.0  Bach  =10.0  Be;  =  0.108407  Sana  :  1.40 


I  :  0.200  4  1.812  T‘2  4  -2.000  7*4  4  0.673  T‘6 


Cfta?  : 

0.001973 

CL  = 

0.043690 

CD  = 

0.009756 

Sl/Sp  : 

2.263438 

V‘(2/3)/Sp  *- 

0.211051 

( L/D)t1b  = 

4.478099 

ib/ibl  : 

0.756626 

7/?i  : 

0.655907 

Sp/spi  : 

0.811647 

b/li  :  0.372398  SSD  = 


0.300503 


Figure  B.4 


Table  B.4 


■ 

Delta  :  8.0 

■ 

8ey  :  0.191*07 

Gaua  :  1.40 

I  :  0.200  ♦  2.457  T‘2  ♦  -3.000  !‘4  ♦  1.061  !*8 

Cfta?  :  0.002009 

CL  =  0.043637 

CD  =  0.009994 

Si/Sp  =  2.399271 

7*(2/3)/S-  :  0.212883 

(L/D)vis  •  4.366130 

Ab/ibi  :  0.S76559 

7/71  =  0.579832 

Sp/Spi  :  0.741174 

b/U  :  0.372398 

SSD  =  0.300503 

_ _ Table  B.5  _ _ 

Phi  :60.0  Delta  :  8.0  Stcb  =10.0  itj  •  C. 101+07  Gam  :  1.40 

I  :  0.300  ♦  0.297  7*2  ♦  C.OOO  T*4  ♦  -0.013  T*« 


CftlT  - 

0.001918 

CL  •- 

0.043599 

CD  : 

0.009593 

S»/Sp  : 

2.107515 

r(2/:j/Sp  j 

0.207946 

(L/D)t1«  : 

4.545040 

Ab/'ibi  = 

0.878447 

f/fi  • 

0.743837 

■  ■  ■  'i 

s? 

0  395?32 

b/ll  = 


C. 409313 


BSD 


0.290238 
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Figure  B. 6 


Table  B.6 


■ 

Phi  =60.0 

_ 

Hfl 

— 

Bach  :10 . 0 

Key  :  0.108+07 

Qsill  :  1.40 

I  s  0.300  ♦  0.411  T*2  ♦  0,000  !*4  ♦  -0.082  8*8 


CftaT  : 

0.001918 

CL  : 

0.043378 

CD  - 

0.009316 

f*(2/3)/Sp  : 

0.207236 

(l/D)?ia  = 

4.379339 

Ab/Abl  s 

0.843229 

f/!l  •- 

0.710444 

Sp/Spi  : 

0.871710 

b/ll  : 

0.409313 

SOD  :  0,290236 

0,290236 


Figure  B.7 


Table  B.7 


Phi  :60.0 

Delta  :  8. 

1 

1 

0.300  ♦  0 

1*2  ♦  -1.000  H  ♦  0.376  1*6 


Cfti?  :  0.001932  CL  :  0.043448 


Sv/Sp  :  2.102324  ?'(2/3)/?p  :  0.203256 


.798890  T/fi  --  0.647930 


0.250238 


b/li  : 


CD  s 


(l/D)tis 


Sp/Spl  : 


0.009450 


4.597888 


.835926 


Figure  B.8 


Table  B.8 


■ 

m 

B 

Bey  =  0.101407 

Sans  :  1.40 

1  S  0.300  4  1.058  T*2  4  .1.000  1*4  4  0.306  T*8 

CL  :  0.043437 

CO  :  0.009431 

St/Sp  -  2.121517 

?‘(2/3)/Sf  :  0.203593 

(L/DM#  :  4,605786 

Ab/A'ol  :  0.783162 

f/Ti  :  0.818193 

8p/Spi  :  0.808807 

b/i*  ;  0.400313 

88D  :  0,200238 
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Figure  B.  9 


b/l»  : 


.409313 


0.290239 


Figure  B.  10 


Table  B.10 


flach  --10.0 

Bey  :  0.101407 

Gain  =  1.40 

1.180  1*2  4  1 

J.OOO  7*4  4  0.011  !‘6 

CftaT  :  0.001891  CL  =  0.043238 


Si/Sp  :  2.077541  ?‘(2/3)/Sp  : 


0.808447  Sin  :  0.821062 


b/lf  •  0.454352  880  :  0.278033 


CD  :  0.009377 


<l/D)'U  =  4.811255 


Sp/Spl  :  0.830038 


Figure  B. 12 
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Figure  B.13 


Table  B.13 


Figure  B.14 


Table  B.14 


■ 

m 

m 

R H  i  0.101+07 

- -  n 

Sam  :  1.40 

I  :  0.400  ♦  2.220  !'2  +  -3.000  T'4  ♦  1.110  T'8 

ClUv  *  0.002045 

Cl  ;  0.042875 

CO  :  0.009654 

S«/Sp  :  2.211403 

f'(2/3)/Sp:  0.197759 

(L/DMi  *-  4.441051 

ib/ibi  *-  0.333009 

f/fi  :  0.371700 

$p/Spi  :  0.593176 

b/li  :  0.454352 

Hgmigg 

APPENDIX  C 


Presented  here  are  configurations  that  have  the  following 

parameters  held  fixed: 

<$j  =  50,  S  *  8,  M  -  10,  Rey  =  1X106 

The  following  parameters  are  varied: 

R  from  0.3  to  0.7 
o 

bj  from  0.0  to  -3.0 

The  values  of  b2  and  b4  are  such  that  equation  2-13  and 
condition  (iii)  holds. 
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Figure  C.l 


Table  C.l 
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Figure  C.2 


Figure  C.3 


Figure  C.4 


Table  C.4 


■ 

Delta  :  8.0 

Hacb  =10.0 

Bey  -  0.101+07 

Gaiia  :  1.40 

X:  0.300  ♦  1.2011*2  ♦  -1.000  n  +  0.339  n 

Cftit  I  0.001907 

CL  :  0.043525 

CD  *•  0.009328 

Sf/Sps  2.170346 

f*(2/3)/Sp  0.209224 

(t/D)»ii  =  4.686131 

Ab/Abl  s  0.027035 

7/fl  s  0.691601 

Sp/Spi  s  0.849130 

b/li  =  0.382059 

SSD  :  0,290238 

Figure  C.5 


Table  C.S 


■ 

Bach  =10.0 

Bar  -  0 . 101+07 

Gana  :  1,40 

I  :  0.300  ♦  1.302  !'2  ♦  -1.200  Y*4  ♦  0.438  T’8 


Cftai  : 

CL  : 

0.043807 

CD  : 

0.009325 

Sl/Sj  : 

2.174980 

V(2/3)/Sp  s 

0.208944 

(l/D)»i«  : 

4.865899 

Ab/Abi  : 

0.818803 

V/?l  = 

0.879508 

Sp/Spi  : 

0.840323 

b/ll  = 

0.382089 

S3D  : 

0.290238 

Figure  C.6 


Table  C.6 


■ 

■ 

Hach  =10.0 

Bey  t  0.10*407 

Gam  :  1.40 

It  0.300  4  1.503  T*2  4  -1.600  T*4  4  0.637  1*6 

CftJT  :  0.001911 

CL  t  0.043470 

CD  ■-  0.009329 

S»/Sp  :  2.187994 

7*(2/3)/Sp  :  0.208499 

(l/D)m  :  4.659933 

ib/ibi  t  0.795958 

7/71  :  0.656134 

Sp/Spi  :  0.822696 

b/l»  :  0.382059 

SSD  t  0.290236 

Figure  C.7 


Table  C.7 


Phi  -50.0 

■ 

Bach  =10.0 

Rejr  *-  0.101407 

Sana  =  1.40 

I  -  0.300  4  1.704  Y*2  4  -2,000  1*4  4  0.838  1*6 

Cftar  :  0.001915 

CL  :  0.043437 

CD  s  0.009345 

Si/Sp  :  2.205945 

?*(2/3)/Sp  =  0.208223 

(l/D)fi«  :  4.548007 

Ab/ibi  :  0.775313 

1/fi  =  0.633565 

Sp/Spi  :  0.804787 

b/lt  =  0.362059 

SSD  s  0.290236 
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Figure  C. 8 


Figure  C.  9 


Table  C.9 


■ 

■ 

Bach  -10.0 

Bey  :  0.101+07 

Oiua  :  1.40 

I  :  0.300  ♦  2.013  T‘2  +  -2.400  Y‘4  +  0.928  Vi 

CftiT  -  0.001926 

Cl  :  0.043415 

CD  -  0.009424 

Sw/Sp  :  2.264482 

7‘(2/3)/Sp  --  0.209875 

(L/D)m  :  4.606875 

ib/ibi  :  0.726901 

7/Vi  =  0.591549 

Sp/Spi  *-  0.762748 

b/l»  :  0.362059 

SSD  :  0.290236 

Table  C.10 


■ 

Delta  :  8.0 

Hack  -10.0 

Bej  s  0 . 108+07 

Sana  =  1.40 

I  :  0.300  ♦  2.312  ?‘2  ♦  -3.000  !*4  ♦  1.229  Y'6 

B2BB1 

Cl  :  0.043382 

— 

CD  :  0.009498 

?'(2/3)/Sp  :  0.210525 

?/fi  =  0,562402 


St/S?  =  2,309484 
Ib/ibi  :  0.696582 
b/l»  =  0.362059 


SSD 


0.290236 


(L/D)Tis  :  4.567349 

Sp/Spi  r  0,735208 


Figure  C.ll 


Figure  C.12 


Table  C.12 


Phi  =50.0  Delta  =  8.0  Mach  =10.0  Ref  =  0. 101+07  Qaaia  :  1.40 


1  =  0.400  ♦  0.642  T'2  ♦  -0.400  Y*4  ♦  0.197  H 


Cf tar  = 

0.001887 

CL  = 

0.043220 

Sw/Sp  - 

2.099222 

7*(2/3)/Sp  = 

0.201318 

ib/ibi  = 

9.809208 

7/71  = 

0.624954 

b/U  = 

9.401898 

SSD  = 

0.276033 

CD  t  0.009145 


(L/D)ria  =  4.726161 


SF/Spi  =  0.824831 


Figure  C. 13 


Table  C.13 


Phi  =50.0 

Delta  =  8.0 

Bach  =10.0 

Uj  -  0. 101*07 

Gaua  =  1.40 

I  =  0.400  ♦  0.750  !*2  ♦  -0.400  T'4  ♦.  0.000  T*6 

Cftar  =  0.001880 

EBMfl 

CD  =  0.009122 

Si/Sp  2.106453 

V ( 2/3) /Sp  =  0.201476 

(L/D)vis  =  4.735709 

Ab/Abi  :  0.781440 

?/H  :  0.600350 

Sp/Spi  --  0.802415 

b/li  :  0.401898 

SSD  :  0.276033 

Figure  C.14 


Figure  C.15 


Table  C.15 


B 

Delta  :  8.0 

■ 

hi  -  0. 101*07 

Gaiaa  :  1.40 

Is  0.400  ♦  1.052  T'2  *  -1.000  7*4  ♦  0.388  7*6 

MllflllllM 

EBHSB 

CD  s  0.009111 

Si/Sp  :  2.113615 

7*(2/3)/Sp  S  0.200432 

(L/Djfis  s  4.734126 

Ab/Abi  -  0.750343 

7/Vi  :  0.566081 

Sp/Spi  :  0.775592 

b/U  s  0.401898 

SSD  :  0.276033 

Figure  C. 16 


Table  C. 16 


Figure  C. 17 


Table  C.17 


■ 

■ 

Each  =10.0 

ley  =  0.101+07 

Sam  i  1.40 

I  =  0.400  ♦  1.562  !‘2  ♦  -1,800  !*4  ♦  0.679  T‘6 

Cftat  :  0.001919 

CL  :  0.043054 

CD  *-  0.009171 

Si/Sp  :  2.167452 

7*(2/3)/Sp  :  0,201086 

(L/D)»ii  :  4.694574 

&b/Abl  --  0.681286 

7/71  =  0.503584 

Sp/Spl  -  0.715070 

b/li  =  0.401898 

SSD  =  0.276033 

Figure  C.18 


Figur*  C. 19 


T«bl«  C.19 


Figure  C.20 


Figure  C. 21 


Bach  -10.0 

ley  s  1 

1.101+07 

Sana  :  1.40 

.500  ♦  0.282  T'2  ♦  0.000  T'4  ♦  0.04?  T*8 


CitjT  : 


Sf/Sp  = 


Ab/Abi  : 


b/lff  r 


Cl  = 

0.042867 

»*(2/3)/Sp  - 

0.193480 

f/fl  S 

0.529250 

SSD  = 

0.258054 

CD  : 


(L/D)t1i  •  4.769799 


Sp/Spi  :  0.768223 
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Table  C.22 


■ 

■ 

Bach  -10.0 

ley  :  0.101+07 

Sana  -  1.40 

X  =  0.500  ♦  0.801  n  ♦  -0.400  T'4  ♦  0.139  !*8 

Cftar  S  0.001877 

CL  :  0.042788 

CD  :  0.008942 

Sa/Sp  :  2.067719 

!'{2/3)/Sp  :  0.192469 

(L/D)»ia  :  4.784766 

Ab/Abl  s  0.704728 

1/n  -  0.483839 

Sp/Spi  :  0.727426 

b/li  s  0.451589 

SSD  s  0.258054 

Figure  C.23 


Table  C.23 


Phi  =50.0 

Delta  :  8.0 

Bach  =10.0 

8er  :  0. 10E+07 

Sana  =  1.40 

I  :  0.500  ♦  0.903  T'2  ♦  -1.000  T‘4  ♦  0.437  1*6 

CftjT  :  0.001094 

Cl  --  0.042716 

| 

Si/Sp  2.075381 

V*(2/3)/Sp  =  0.191391 

(L/D)vis  s  4.773218 

Ab/Abi  :  0.673631 

T/Ti  =  0.452648 

8p/Spi  :  0.699733 

b/li  :  0.451589 

8SD  s  0.258054 

Figure  C. 24 

Table  C.24 

Phi  =50.0  Delta  *  8.0  Sach  =10.0  ley  =  0. 181*07  Gaaia  =  1.40 

1  :  0.500  ♦  1.212  n  *  -1.400  1*4  *  0.529  T*6 

Cftar  *.  0.001915  CL  :  0.042884  CD  =  0.008987 

S»/Sp  =  2.105814  f‘(2/3)/Sp  *.  0.192181  (L/D)t1»  x  4,747234 

kW  ■  0.825218  f/fi  :  0.414169  Sp/8pi  =  0.656853 

b/lf  =  0,451589  SSD  s 


.256054 


Figure  C. 25 


Table  C.25 


Figure  C.26 


Table  C.26 


Figure  C.27 


Table  C.27 


Phi  :50.0 

EB 

m 

ley  :  0.101+07 

6aiaa  =  1.40 

I  :  0.500  ♦  1.815  1*2  ♦  -2.600  T‘4  ♦  1.126  T*6 

Cftar  :  0.001955 

Cl  :  0.042571 

CD  :  0.009131 

S>/Sp  =  2.172482 

r(2/3)/Sp  :  0.193128 

(l/D)ris  :  4.662208 

Ab/Abi  :  0.563284 

f/fi  :  0.363326 

Sp/Spl  :  0.598916 

b/lw  :  0.451589 

SSD  :  0.258054 

Figure  C.28 


Bach  =10.0 

fiey  :  ( 

). 106+07 

Gaaaa  =  1.40 

Is  0.500  ♦  2.016  Y‘2  ♦  -3.000  ?*4  ♦  1.325  T*6 


CO  :  0.009198 


(L/D)vis  :  4.626232 


Sp/Spi  s  0.579266 


Cftar  = 

0.001969 

CL  s 

0.042552 

St/Sp  = 

2.201684 

9’(2/3)/Sp  = 

0.194077 

ib/Abl  = 

0.542639 

9/71  = 

0.348142 

b/ll  s 

0.451589 

SSD  : 

0.258054 

Figure  C.29 


Table  C.29 


Phi  =50.0  Delta  :  8.0  Bach  =10.0  Bey  :  0 . 108+07  Sana  :  1.40 


X  :  0.600  +  0.143  1*2  ♦  0.000  7*4  ♦  0.096  7*6 


CD  --  0.008833 


(L/D)Tia  =  4.806467 


Sp/Spi  :  0.692483 


Cftar  = 

0.001841 

CL  = 

0.042447 

S»/Sp  : 

2.06106? 

V(2/3)/Sp  *- 

0.184276 

ib/Ibl  = 

0.676428 

7/71  = 

0.421012 

b/li  : 

0.616302 

c n 

»  1 

0.236489 

Figure  C.30 


Table  C.30 


■ 

B 

Hacb  -10.0 

Her  s  0.101+07 

Gaiaa  :  1.40 

Is  0.600  +  0.251  T‘2  ♦  0.000  T‘4  t  -0.0111*6 

Cftat  s  0.001851 

CL  s  0.042412 

CD  ;  0.008800 

St/Sp  :  2.045426 

Y*!2/3)/Sp  :  0.184068 

(L/D)?ia  :  4.819488 

Ab/Abi  :  0.648661 

y/!i  :  0.399285 

Sp/Spi  s  0.669200 

b/l«  s  0.515302 

SSD  s  0.236489 

Figure  C.31 


j 

i 


Table  C.31 


1 

■ 

1 

Bach  =10.0 

Be i  •  0,101+07 

Gam  :  1.40 

I  :  0.600  ♦  0.452  !*2  ♦  -0.400  T‘4  ‘  0.188  1*6 


Cftav  : 

0.001867 

CL  : 

0.042357 

CD  : 

0.008801 

Sfl/Sp  : 

2.043804 

T*(2/3)/Sp  s 

0.183058 

(L/D)tU  S 

4.812601 

Ab/Abi  : 

0.628016 

7/71  = 

0.379612 

Sp/Spl  : 

0.650607 

b/li  -- 

0.515302 

m 

0.236489 

Figure  C.32 


Table  C.32 


Bach  -10,0 


S«7  2  0. 101*07  Gun  :  1.40 


Phi  :S0.l  Delta  :  1,0 


I*  0,000  ♦  1.080  1*2  ♦  -1.400  !‘4  ♦  0.5711*0 


CftlT  s  1 

1.001020 

Ot/Op  *  1 

1.075410 

Ab/ibi  :  I 

1.540508 

b/li  i  1 

1.515002 

0.042220 


0.102000 


0.017300 


0.200400 


CD  :  0.008071 


(L/D)»i»  s  4.758073 


8p/8pi  :  0.570470 
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Figure  C.34 


Table  C.3A 


Phi  =50.0  Delta  =  8.9  Hich  =10.0  9ay  =  0.101407  Quit  :  1.40 


I  :  0.800  i  1.284  T*2  4  -1.800  T‘4  4  0.777  !*i 


Cftl? 


Sa/Sp  : 


ib/M  : 


b/li  -* 


CL  : 

0.042183 

T*(2/3)/5p  : 

0.162802 

?/fl  = 

0.301864 

S8D  = 

0.236489 

CD  :  0.008919 


(L/3)rla  :  4.729744 


Sp/Spt  :  0.558982 


Figure  C.35 


Figure  C.36 


Table  C.36 


Phi  =50.0  Delta  =  8.0  Mach  =10.0  hj  :  0.101+07  Gam  =  1.40 


I  :  0.800  ♦  1.868  T'2  ♦  -2.800  T*4  ♦  1.175  1*8 


Cftar  : 

0.001980 

Cl  = 

0.042127 

Sl/Sp  : 

2.138979 

7*(2/3)/Sp  = 

0.184060 

Ab/Abi  s 

0.486572 

f/fi  = 

0.273058 

b/la  = 

0.515302 

SSD  - 

0.236489 

CD  :  0.009043 


(l/D)fie  =  4.658574 


Sp/Spi  :  0.519468 


Figure  C. 37 


Figure  C.38 


Figure  C.  39 


Table  C.39 


■ 

■ 

Bey  s  0.101+07 

Qtua  -  1.40 

I  :  0.700 

♦  0.411  T‘2  ♦  -0.400  T*4  ♦  0.130  T*6 

Cftar  = 

0.001075 

CL  :  0.041883 

CD  : 

0.008686 

S»/8p  = 

2.033808 

?*(2/3)/Sp  =  0.173457 

(l/D)?i»  : 

4.821683 

Ab/Abi  : 

0.52353$ 

T/Ti  :  0.270647 

Sp/Spi  : 

0.547973 

b/ll  = 

0.599945 

SSD  -  0.211550 

Figure  C. AO 


Table  C.AO 


■ 

■ 

Bach  =10.0 

Bey  :  0.101+07 

Gaaia  =  1.40 

I  :  0.700  ♦  0.712  T*2  ♦  -1.000  T‘4  ♦  0.429  T‘6 

Cftar  :  0.001911 

CL  =  0.041807 

CD  --  0.008740 

Sf/Sp  :  2.044167 

r<2/3)/Sp  :  0.172623 

(l/D)»ia  :  4.783686 

4b/Abi  -  0.492698 

?/fi  s  0.247583 

Sp/Spi  :  0.518875 

b/lf  s  0.599945 

SSD  =  0.211550 

Figure  C.41 


Table  C.41 


B 

m 

Mach  :10.0 

hj  :  0.108*07 

6sus  :  1.40 

1  =  0.700  +  0,913  7*2  *  -1.400  7*4  +  0.828  7*6 

CL  :  0.041763 

CD  s  0.008788 

Si/Sp  :  2.056345 

7*(2/3)/Sp  :  0.172440 

(L/D)tU  --  4.752065 

ib/ibi  s  0.472053 

7/71  :  0.233263 

Sp/Spi  :  0.499199 

SSD  : 


b/lf  : 


0.599945 


0.211550 


Figure  C.42 


Figure  C.43 


Table  C.43 

Pbi  =50.0  Delta  =  8.Q  Bach  =10.0  Bey  :  0 . 108+07  Gaua  :  1.40 


I  :  0,700  4  1.315  7*2  4  -2.200  7‘4  4  1.026  T*6 


Cftar  = 

0.001987 

CL  = 

0.041693 

CD  : 

0.008916 

Sl/Sp  : 

2.091982 

y*(2/3)/Sp  -- 

0.173211 

(l/D)n»  = 

4.676222 

Ab/Abi  = 

0.430764 

m  -- 

0.207317 

Sp/Spi  = 

0.459405 

SSD  : 


b/li  = 


0.599945 


0.211550 


125 


Figure  C.44 


Table  C.44 


Figure  C.45 


APPENDIX  D 


Presented  here  are  conf igurations  that  have  the  following 
parameters  held  fixed! 

$,  -  40,  6-8,  M  -  10,  Rey  -  IX 10* 

The  following  parameters  are  varied! 

Rfl  from  0.4  to  0.8 
bj  from  0.0  to  -2.0 

The  values  of  bj  and  are  such  that  equation  2-13  and 
condition  (iii)  holds. 
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Figure  D.l 


Table  D.l 


hi  *40.0 

Delta  *  0.6 

laeb  *10,0 

Iff  *  0.101*07 

Oiiii  *  1.40 

I  *  0.400  ♦  1.070  7*2  4  0,000  T'4  ♦  -0,420  T’O 

Cftiv  *  0.001072 

Cl  *  0.04)200 

CD  *  0.000160 

It/Op  *  2.220440 

?*(2/2)/0p  *  0.210020 

<l/D)»U  *  4.720020 

mama 

7/71  *  0.040400 

Op/Spl  *  0,020101 

b/li  *  0,217222 

000  *  0.27001) 

J> I 


Figure  D.2 


T»bl*  D.2 


Figure  D.3 


Table  D.3 


m 

m 

Bach  =10.0 

Bey  :  0.101+07 

Oaaaa  =  1.40 

I  :  0.400  ♦  1.778  7*2  ♦  *2.000  7*4  ♦  0.977  7*8 

Cfttv  :  0.001879 

CD  :  0.009145 

Si/Sp  :  2.240033 

7*(2/3)/Sp  =  0.209406 

(L/D)t1s  =  4.719671 

ib/Abi  :  0.764798 

7/71  :  0.585486 

Sp/Spi  :  0.779723 

b/l»  :  0.337233 

SSD  :  0.278033 

Figure  D.4 


Table  D.4 


Figure  D. 5 


Table  D.5 


Phi  =40.0  Delta  =  8.0  Mach  =10.0  ley  =  0,101+07  Sim  =  1.40 

I  :  0.500  +  0.858  T‘2  ♦  0.000  T*4  ♦  *0.205  T‘8 

Cftar "=  0.001859  |  CL  =  0.042861  CD  :  0.008015 

!*(2/3)/Sp  :  0.201520  (L/D)ri«  =  4.607703 

T/Tl  =  0.521240  Sp/Spi  =  0.740606 


Si/Sp  :  2.151700 
Ib/lhl  ;  0.738083 
b/lt  :  0.378020 


0.258054 


Figure  D.  6 


Table  D.6 

Phi  MO.O  Delta  :  8.0  list  =10.0  ley  :  0.101*07  Oiua  :  1.40 

I  •-  0.500  ♦  1.212  !*2  ♦  *1.000  !*4  *  0.411  T*0 

Cftiv  :  0.00110$  CL  *  0.042m  CD :  0.001012 

Sv/Sp  :  2.154411  r(2/3)/3p  :  0.200740  (L/D)vli  :  4.101750 

ib/Abl  :  0.711191  f/U  :  0.49211)  Sp/Spl  *  0.725149 

b/li  S  0,178020  1SD 


.251054 


Figure  D.7 
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Figure  D.8 


Table  D.8 


Figure  D.  10 


Table  D.10 


■ 

m 

m 

Bey  :  0.10B+07 

Gaaaa  =  1.40 

I  S  0.600  ♦  1.45?  T‘2  ♦  -2.000  1*4  ♦  1.053  T*6 

Cf tar  :  0.001884 

CL  :  0.042278 

CD  :  0.008783 

Ss/Sp  :  2.127193 

7"(2/3)/Sp  x  0.191596 

(L/D)»is  i  4.813523 

Ab/ibi  =  0.584419 

7/71  :  0.347349 

Sp/Spi  x  0.601709 

b/li  s  0.432390 

SSD  :  0.236489 
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Figure  D.ll 


Table  D.ll 


Figure  D. 12 


Table  D.12 


n 

■ 

Bach  =10.0 

Bey  =  0.108407 

Gaua  =  1.40 

I  =  0.700  4  0.891  n  4  -1.000  Y‘4  4  0.487  Y‘S 

| 

CL  :  0.041880 

CO  =  0.008588 

Str/Sp  :  2.085479 

?'(2/3)/Sp  =  0.181381 

(L/0)fl»  :  4.878789 

Ab/Abl  :  0.530814 

7/71  •  0.276388 

Sp/Spi  :  0.545752 

b/li  :  0.503414 

SSD  =  0.211550 

Figure  D. 13 


Table  D.13 


Figure  D. 14 


Table  D. 14 


■ 

m 

m 

Bey  :  0.101407 

Gian  *  1.49 

I  :  0,800  ♦  9.222  T*2  ♦  0.099  T*4  ♦  9.120  T*0 

Cftar  :  0.091014 

CL  :  9.041510 

CD  *  9.000492 

Sv/Sp  :  2.931092 

!*(2/3)/9p  :  0.179501 

(l/D)fli  :  4.941123 

Ob/Abl  •-  0.490075 

7/Ti  :  9.225900 

Sp/Spi  :  0.507002 

aana 

OSS  :  0.100409 

Figure  D. 15 


APPENDIX  E 


Presented  here  are  configurations  that  have  the  -following 

parameters  held  fixed: 

*=  30,  S  =  S,  M  =  10,  Rey  =  1X106 

The  following  parameters  are  varied: 

R„  from  0.5  to  0.9 
o 

b|  from  0.0  to  -15.0 

The  values  of  bj  and  b^  are  such  that  equation  2-13  and 
condition  (iii)  holds. 
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Figure  E.l 


Table  E.l 


Phi  =30.0  Delta  =  8.0  Bach  =10.0  Her  =  0.101+07  Gain  =  1.40 


I  :  0.500  ♦  2.010  1*2  ♦  0.000  T’4  ♦  -2.943  T*« 


C 


Sl/Sp  : 


ib/Abi  : 


b/lf  2 


CL  : 

0.042870 

l*(2/3)/Sp  2 

0.217357 

f/fi  2 

0.515734 

SSD  2 

0.258054 

0.009194 


(L/D)fie  :  4.667721 


Sp/Spi  :  0.732560 


Figure  £.2 


T4bl«  E.3 


Phi  =30.0 

Pilti  s  8.9 

lick  *10.0 

hr  s  0.101+07 

61m  *  1.40 

is  o.ioo  ♦  i.ieon  ♦  *10.000 ri  ♦  s,?66n 

Cftiv  s  0.001061 

Cl  s  0. 942671 

CP  :  0.000692 

Si /It  :  2.600001 

?*(2/l)/$p  :  0 .219706 

(l/DM*  :  4.646412 

ib/Obl  :  0.620006 

f/Tl  :  0.410616 

Sp/Spl  :  0.634407 

b/li  s  0.294764 

SIP  *  0.261064 

Figure  E. A 


Figure  E.5 


Bach  -10.0 


Rei  :  0.101+07 


Gaiaa  =  1.40 


Figure  E.7 


Figure  E.8 


Table  E.8 


Figure  E.9 


Table  E.9 


Pbl  =30.0  Delta  :  8.0  Bach  =10.0  8 tj  :  0.108+07  Gana  :  1.40 


0.700  ♦  1.318  1*2  ♦  0.000  2*4  ♦  -1.667  7*6 


Cftll  : 

0.001840 

CL  s 

0.041957 

CD  = 

0.008657 

S»/Sp  : 

2.183216 

n 

cv 

to 

<■0 

CM 

0.196271 

iL/D)fie  = 

4.846644 

Ab/Abi  s 

0.858292 

T/fl  : 

0.300542 

Sp/Spi  = 

0.565997 

b/lf  : 

0.391586 

SSD  : 

0.211550 

Figure  E.10 


Table  E.10 


Phi  =30.0 

Delta  -8.0 

Bach  =10.0 

hi  :  0.101+07 

Giua  :  1.40 

I  :  0.700  ♦  2.389  1*2  ♦  -5.000  1*4  ♦  4.187  T*6 

Cftar  :  0.001871 

CL  =  0.041834 

CD  :  0.008763 

Si/Sp  :  2.221092 

V* (2/3)/Sp  =  0.196094 

(l/D)fie  =  4.774157 

ib/ibi  :  0.507235 

7/71  =  0.261449 

Sp/Spi  :  0.516250 

b/li  :  0,391586 

SSD  =  0.211550 

Figure  E. 11 


Figure  E. 12 


Table  E. 12 


■ 

Delta  :  0.0 

Bach  slD.O 

Bey  ?  0. 101*07 

Sana  =  1.40 

I  *.  0.700  ♦  4.512  T*2  ♦  -15.000  T'4  ♦  15.035  T'6 

Cftir  :  0.001036 

CL  :  0.041715 

CD  :  0.000240 

Si/Sp  s  2.424003 

7*(2/3)/6p  :  0.203740 

(L/D)t1(  :  4.510335 

ib/ibi  :  0.405110 

7/T1  :  0.200114 

8p/8pi  :  0.415750 

b/lf  =  0.301566 

■At  til 


Figure  E. 13 


Figure  E.1A 


Table  E.U 

Phi  =30.0  Delta  =  8.0  Mach  =10.0  Her  =  0.101+07  6aua  =  1.40 

I  =  0.800  t  2.038  1*2  ♦  -5.000  T*4  ♦  4.804  T'6 

Cftar  =  0.001883  CL  =  0.041357  CD  =  0.008578 

St/Sp:  2.140818  7‘(2/3)/$p  =  0.183749  (l/D)fi»  =  4.822828 

ib/4bi  =  0.423858  7/71  =  0.181456  Sp/Spi  =  0.431872 

b/lt  =  0.488550  S8D  =  0.183480 


Figure  G. IS 


Figure  E.16 


Figure  E.17 


Table  E.17  _ 


Figure  E.  18 


Table  E.18 

Phi  =30.0  Delta  ;  8.0  Hack  =10.0  Bey  :  0.101+07  Sana  :  1.40 

X  :  0.800  +  1.888  T'2  +  *8.000  T*4  +  5.442  T'8 

Cftav  :  0.001898  CL  :  0.040875  CD  =  0.008433 

Se/Sp  =  2.079087  ?‘(2/3)/Sp  :  0.188835  (L/D)ila  :  4.847222 

ib/ibi  :  0.340499  7/71  :  0.116021  Sp/Spi  =  0.347197 

b/li  =  0.583168  88D  =  0.152450 


Figure  E. 19 


Table  E.19 


Phi  =30.0 

■ 

8ich  =10.0 

Bey  s  0.108+07 

Oaaia  =  1.40 

I  :  0.900  ♦  2.759  T‘2  ♦  -10.000  T*4  ♦  11.280  1*8 

|| 

CL  *-  0.040799 

CD  s  0.008725 

Si/Sp  S  2.156927 

7*(2/3)/Sp  :  0.173278 

(1/D)t1s  --  4.676111 

Ab/4bi  :  0.289528 

T/!i  =  0.094822 

Sp/Spi  =  0.296699 

b/lt  -  0.583188 

SSD  -  0.152450 

APPENDIX  F 


Presented  here  are  configurations  that  have  the  -following 

parameters  held  fixed: 

=  20,  8=8,  M  =  10,  Rey  =  1X106 

The  following  parameters  are  varied: 

R  f rom  1.0  to  1.1 
o 

b^  from  8.0  to  -8.0 

The  values  of  b2  and  b^  are  such  that  equation  2-13  and 
condition  (iii)  holds. 
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Figure  F.l 


Table  F • 1 


Figure  F.2 


Table  F.2 
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Figure  1’ .4 


Figure  F.5 


Table  F.5 


■ 

m 

Bach  =10.0 

I«7  •  0 . 101+07 

Suit  i  1.40 

I  :  MOO  ♦  0 .777  T'2  ♦  0,000  T*4  ♦  -3.210  T‘0 


Cftav  : 

0,001016 

CL  ■ 

0.040021 

CD  * 

0.0(0021 

J»/8p  • 

2.030001 

rm/b  * 

0.103377 

(l/D)»i»  : 

4.9(0471 

Ib/iftl  : 

0.201037 

f/fl  : 

0.037003 

tlM 

0.203130 

BB 

0.701022 

SOD  J 

0.002074 

Pigur*  P.6 


T«bl«  P.6 


APPENDIX  6 


Presented  here  are  configurations  that  have  the  following 
parameters  held  fixed: 

-  10,  $  -  8,  M  «  10,  Rey  «  1X106 

The  following  parameters  are  varied: 

Rq  f  rom  1,0  to  1,1 
from  120.0  to  -120.0 

The  values  of  bj  and  are  such  that  equation  2-13  and 
condition  (iii)  holds. 
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Figure  G.l 


Teble  G.l 

Phi  =10.0  Delta  =  8.0  Itch  =10.0  lay  =  0,101*07  G«».  ;  1.40 

1  :  1.000  *  5. 088  T*2  *  120,000  1*4  *  *2138,0  1*8 

Cftir  :  0.001814  Cl  **  0.040848  CO  =  0.008108 

8»/8p  =  2.828714  ?*(2/8)/8p  =  0,222710  (l/D)»U  :  4.482078 


ib/lhl 

: 

.282392 

f/fi  : 

,083182 

8p/8pi  : 

1.293198 

b/ll  = 

.288121 

880  :  0 

,118747 

Figure  G.2 


Table  G.2 


Phi  *.10.0 

Delta  :  0. 

I 

Bach  -10 . 

0 

Oej  :  0.101407 

Still  *  1.40 

I 

1.000  4 

8.2oo  n  ♦ 

1.000  T*4  4  -077.0  T'O 

Cftir  :  0,0913)1  CL  :  0.040500  CO  »  0.009214 
S»/8p  :  2.0(34)4  f‘(2/3)/0p  :  0.222003  (L/0)vli  :  4.305521 
ib/Abl  *  0.274031  ?/?i  :  9.075400  Sp/lpi  s  0.274054 
b/li  :  0.200121  000  :  0.110747 


Figure  G.3 


Table  G.3 


m 

Hicb  MO.O 

lif  :  0.101*07 

Oiiii  *  1,40 

I  :  1,000  ♦  11.300  !*2  W20.000  T*4  ♦  101,0  T*0 


CftlT  : 

0,001844 

mm 

0.040402 

CP : 

0.000307 

Sl/Sp  : 

2,731030 

»*(2/J)/0p  * 

0.224272 

(L/D)fll  : 

4.310472 

Ab/ibl  : 

0.200343 

f/fl  : 

0.000000 

mm 

m 

b/l»  J 

0,200121 

000  : 

0.110747 

Figure  G.A 


Table  G.A 


Figure  G.S 


Table  G.5 


■ 

a 

a 

Bey  :  0. 101*07 

inn  •  1.40 

X:  UOO  ♦  5.295  T*2  ♦  0.000  1*4  ♦  -014.0  1*0 

Cltii  :  0.001041 

CL  :  0.040009 

CO  :  0.000000 

it /it  2.3210(9 

f*(2/3)/0p  :  0.190000 

(L/D)rli  :  4.(7(000 

Ib/lbi  :  0.107100 

7/fl  :  0.034(00 

Sp/Spi  :  0.107090 

b/li  :  0.3900)0 

Figure  G.6 


APPENDIX  H 


Here  in  appendix  H  we  show  several  different 
conf igurations  in  a  three-dimensional  setting.  The 
configurations  shown  are  repeats  of  earlier  configurations  in 
the  above  appendices. 

These  3~D  views  are  presented  so  that  the  reader  gets  a  better 
perception  of  the  conf igurations  shown  and  even  ones  that  are 
shown  only  in  the  2-D  views  of  earlier  appendices. 
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Figure  H. 1  ( C . 3 ) 
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Figure  H.2  (C.45) 


Figure  H.5  (E.10) 


F igure  H . 8  (  G  .  6  ) 


